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1.0  SUMMARY 


A  theoretical  research  investigation  has  been  undertaken  of  high  energy  chemical  species 
that  are  attractive  candidates  for  advanced  propulsion  systems.  The  tasks  of  this  research 
program  to  date  have  included  the  follovving  studies:  1)  high  symmetry  and  hypervalent 
structures  as  advanced  oxidizers;  2)  azide-like  species  such  as  azidamine  (Ny)  and  nitryl  azide 
(NO2N3);  3)  nitrated  cyclic  and  polycyclic  compounds,  as  further  developments  of  cyclic 
nitramines  such  as  the  cyclotetramethylenetetranitramine  (HMX)  and  polycyclic  nitramine 
(PCN)  family  of  propellants;  4)  an  extensive  analysis  of  dinitramide  compounds  as  potential 
chlorine-free  oxidizers;  5)  nitrated  ring  structures  including  nitrozoles  and  azididines;  and  6) 
substituted  prismanes.  The  significant  findings  in  these  six  areas  of  high  energy  materials  are 
described  in  the  folloviing  paragraphs. 

Studies  of  oxygen  hexafluoride  (OF5)  indicate  that  this  high  symmetry  stracture  may  exist  as 
a  gas  phase  species  exhibiting  mainly  covalent  bonding.  However,  there  is  a  low  transition  state 
barrier  for  formation  of  the  ionic  oxygen  hexafluoride  (F'OF5'^  structure,  which  has  low  energy 
content.  The  chemistry  of  this  compoimd  thus  resembles  that  for  nitrogen  pentafluoride  (NF5). 

It  may  be  of  interest  to  study  oxygen  tetrafluoride  (OF4)  and  other  mixed  oxidizers,  but  the 
potential  for  significant  specific  impulse  (Igp)  improvement  from  this  class  of  oxidizers  appears 
to  be  limited. 

All  of  the  new  azide-like  structures  that  were  examined  exhibit  very  large  positive  heats  of 
formation.  The  chemistry  of  the  diazidamide  anion  parallels  that  for  the  dinitramide  anion  but 
exhibits  much  higher  energy  content.  The  fundamental  problem  with  azidamines  is  their  low 
energy  barrier  for  nitrogen  elimination.  This  class  of  compounds  (azides)  is  chemically 
interesting,  highly  energetic,  but  exhibits  limited  structural  stability,  even  offering  explosive 
character.  Further  research  in  this  area  is  warranted,  however,  since  synthesis  techniques  for 
azide-like  structures  are  often  applicable  to  other  classes  of  compoimds  that  exhibit  low 
decomposition  barriers. 

Nitramine  chemistry  is  presently  a  dynamic  area  for  synthesis  of  new  compounds. 
Calculations  of  hexanitroisowurtzitane  (CL20)  confirm  the  structure  of  this  molecule  and  verify 
the  high  positive  heat  of  formation.  Based  on  these  studies  of  CL20,  several  new  polycyclic 
nitramine  structures  were  analyzed.  A  second  compound  of  this  type  can  be  formed  by  carbon- 
bridging  on  two  HMX  molecules.  Nitrogen-bridging  was  examined  on  two 
cyclotrimethylenetrinitramine  (RDX)  molecules  and  a  stable  pentadecane  stmcture  was  found: 
2,4,6,8,10,12,13,14,15-nonanitra-2,4,6,8,10,12,13,14,15-nonaazatetracyclon  (5.5.1. 
pentadecane.  This  is  the  first  polycyclic  nitramine  that  was  analyzed  that  exhibits  a  positive 
oxygen  balance  (+3  O  atoms).  Synthesis  activities  on  these  new  polycyclic  nitramines  appears  to 
be  warranted. 

Several  new  compounds,  mainly  alkali  and  ammonium  salts,  were  synthesized  that  exhibit 
the  basic  dinitramide  anion  [N(N02)2']  structure.  An  extensive  study  of  this  new  inorganic 
anion  and  several  new  dinitramide  compounds  was  imdertaken.  Several  new  compounds  were 
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examined  for  structural  stability  including  trinitramine  [N(N02)3]  and  hydroxyammonium 
dinitramide  [HADN].  This  dinitramide  group  can  also  be  realized  by  systematic  nitration  of  a 
starting  compoimd  such  as  hydrazine.  Three  vibrationally  stable  compounds:  1,1- 
dinitrohydrazine,  1,2-dinitrohydrazine  and  tetranitrohydrazine  were  found.  This  last  compound 
should  be  a  leading  candidate  for  synthesis  activities  of  advanced  oxidizers. 

Studies  of  nitrated  ring  structures  focused  on  the  structure,  stability  and  energy  content  of 
1,3,3-trinitroazididine  [TNAZ].  Beyond  TNAZ,  calculations  were  undertaken  on  several 
compoimds  that  were  suggested  or  already  S5mthesized  by  SRI,  International  (SRI).  They 
include  l-amino-3-nitro- 1,2,4-triazole;  1 -amino-3, 5-dinitro- 1,2,4-triazole;  l-fluoro-3-nitro- 1,2,4- 
triazole;  and  l-amino-3-nitrotetrazole.  Finally,  l,3-dinitro-l,3-diazacyclobutane,  which  is  the 
smallest  member  in  the  RDX,  HMX  nitramine  series  was  examined.  The  studies  of  TNAZ 
suggested  that  this  compound  should  exhibit  stability.  It  was  found  that  this  compound  has  a 
molar  heat  of  formation  similar  to  RDX  and  surprisingly  little  strain  energy  in  the  ring.  The  name 
HEDMX  is  suggested  as  an  identifier  for  this  new  cyclic  nitramine,  which  is  a  prime  target  for 
experimental  synthesis. 

The  last  series  of  compounds  that  was  studied  were  [n]-prismanes.  Calculations  indicate  that 
both  the  nitroprismanes  and  dinitroprismane  structures  are  stable.  Nitration  lowers  the  heat  of 
formation  by  ~  10%,  but  these  compounds  exhibit  high  heats  of  formation  (+100  to  +125 
kcal/mol)  in  any  case,  owing  to  their  high  strain  energy.  Further  synthesis  and  compound 
characterization  activities  should  be  encouraged  for  this  class  of  strained  hydrocarbons. 
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2.0  INTRODUCTION 


This  report  docxmients  the  work  completed  by  The  United  Technologies  Research  Center 
(UTRC)  on  Contract  F0461 1-90-C-0009,  "Theoretical  Research  Investigation  of  High  Energy 
Chemical  Species,"  for  the  United  States  Air  Force  (USAF)  Operating  Location  (OL-AC) 

Phillips  Laboratory  (PL),  Propulsion  Directorate,  Edwards  AFB,  California.  It  was  intended  that 
the  theoretical  studies  under  this  contract  would  be  closely  coupled  with  experimental  synthesis 
activities  on  many  of  these  compounds  that  were  ongoing  at  SRI  International,  Rockwell  Science 
Center,  Rocketdyne,  and  Phillips  Laboratory.  The  objectives  of  this  research  program  were  to 
identify  those  light  element  molecular  species  that  have  sufficient  internal  energy  to  be  useful  in 
advanced  chemical  propulsion  systems;  to  characterize  their  molecular  structure,  energy  content 
and  structural  or  radiative  stability;  to  examine  the  probable  kinetic  routes  for  decomposition  of 
these  species;  to  characterize  the  vibrational  and  electronic  spectra  and  other  static  molecular 
properties;  and  to  propose  possible  synthesis  routes  for  their  formation  and  subsequent  storage. 

A  goal  of  this  theoretical  investigation  was  to  furnish  fundamental  data  and  to  provide  a 
better  imderstanding  of  the  possibilities  of  forming  chemical  compoxmds  that  represent  advanced 
high  energy  sources.  It  was  also  intended  that  this  theoretical  research  effort  be  carried  out  with 
sufficient  flexibility  to  aid  in  the  analysis  and  guidance  of  any  concurrent  experimental  programs 
of  synthesis  of  high  energy  molecular  species.  To  this  end,  the  tasks  of  this  research  program  to 
date  have  included  the  following  studies:  1)  high  synunetry  and  hypervalent  structures  as 
advanced  oxidizers;  2)  azide-like  species  such  as  N7  and  NO2N3;  3)  nitrated  cyclic  and 
polycyclic  compoimds,  as  further  developments  of  cyclic  nitramines  such  as  the  HMX  and  PCN 
family  of  propellants;  4)  an  extensive  analysis  of  dinitramide  compounds  as  potential  chlorine- 
free  oxidizers;  5)  nitrated  ring  structures  including  nitrozoles  and  azididines;  and  6)  substituted 
prismanes. 

The  research  program  was  formulated  in  three  phases:  The  first  phase  was  devoted  to 
screening  potential  advanced  propellants  in  the  six  groups  described  above.  On  the  basis  of  these 
studies,  several  candidate  species  were  selected  for  more  accurate  calculations  and  predictions  of 
their  thermochemistry.  The  final  phase  involved  studies  of  molecular  properties  of  several 
selected  advanced  propellant  species  and  an  analysis  of  their  potential  energy  reaction  surfaces  in 
regions  defining  the  formation  and  decomposition  of  these  species. 

During  the  course  of  this  research  investigation,  several  new  types  of  compounds  were 
suggested  based  on  the  results  of  the  screening  studies.  In  particular,  exploratory  studies 
indicate  stability  for  trinitramine  [N(N02)3],  and  for  several  salts  formed  from  the  tetranitramide 
cation  [N(N02)4‘^].  In  addition,  preliminary  studies  indicate  stability  for  tetranitrohydrazine 
[^2(1402)4],  which  has  potential  as  a  highly  energetic  oxidizer  material.  Exploratory  studies  also 
found  stability  for  several  [n]-prismane  and  nitro-prismane  compoxmds.  Unsubstituted 
triprismanes,  cubane  and  pentaprismane  structures  have  been  synthesized.  Working  in 
collaboration  with  Dr.  Karl  Christe,  UTRC  identified  and  characterized  a  class  of  azide 
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compoimds  known  as  azidimines  and  carried  out  further  characterization  of  dinitramide 
structures.  An  extensive  set  of  calculations  was  also  carried  out  on  many  of  the  nitrated  triazole 
and  pentazole  compounds  suggested  or  synthesized  by  SRI. 

In  addition  to  the  compounds  described  above,  UTRC's  theoretical  data  on  bonding  and  the 
relative  stability  of  several  new  compounds,  particularly  data  derived  from  various  nitramide 
forms,  had  suggested  a  new  class  of  inorganic  compounds  based  on  the  dinitroboride  anion 
[B(N02)2']-  In  preliminary  studies,  stability  was  found  for  these  structures  which  have  neither 
been  synthesized  or  discussed  in  the  extensive  literature  on  boron  compounds.  Detailed  studies 
of  this  new  class  of  compounds,  however,  led  to  low  barriers  for  rearrangement  to  imenergetic 
forms  or  to  decomposition. 

The  calculated  results  and  pertinent  discussion  are  presented  in  Section  3.  A  concluding 
summary  and  recommendations  based  on  the  results  of  this  program  are  given  in  Section  4. 
Appendix  A  contains  a  list  of  all  technical  presentations  and  publications  that  have  resulted  from 
this  research  program.  Abstracts  of  the  pertinent  published  papers  are  then  collected  in 
Appendix  B. 
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3.0  DISCUSSION  OF  RESULTS 


This  research  program  was  divided  into  studies  of  molecular  species  that  can  be  categorized 
into  several  groups  as  follows: 

1 .  high  symmetry  hypervalent  structures, 

2.  azide-like  structures, 

3 .  nitrated  cyclic  and  polycyclic  compounds, 

4.  dinitramide  compounds, 

5.  nitrated  ring  structures, 

6.  substituted  prismanes. 

Detailed  reports  on  many  of  these  studies  have  been  published  in  the  scientific  literature,  as 
shown  in  Appendix  A.  A  brief  summary  of  the  calculated  results  is  presented  in  Sections  3.1 
through  3.6,  as  categorized  by  the  groups  listed  above  in  this  paragraph. 

3.1  High  Symmetry  Hypervalent  Structures  (OF5) 

Ab  initio  calculations  of  OFg  were  carried  out  at  both  the  SCF  and  MP2  levels  of  theory, 
using  several  different  basis  sets.  Calculations  were  performed  using  both  the  Gaussian  and 
CADPAC  series  of  programs,  the  latter  mainly  for  second  derivative  analysis  at  the  MP2  level 
for  this  closed  shell  structure.  Optimized  OFg  geometries  were  calculated  at  the  SCF  and  MP2 
levels  of  theory.  Harmonic  vibrational  frequencies  were  subsequently  computed  at  the  SCF  and 
MP2  stationary  points. 

Results  of  the  geometry  optimization  calculations  are  shown  in  Table  1 .  All  calculations 
predict  an  octahedral  structure.  Oh.  The  optimum  bond  length  for  OFg  appears  to  be  only  weakly 
dependent  on  the  chosen  basis  set  at  the  HF  level  of  theory.  However,  at  the  MP2  level  of 
theory,  there  is  a  significant  lengthening  of  the  bond  if  diffuse  functions  are  included  in  the  basis 
set.  If  the  optimum  bond  length  becomes  greater  than  ~  1 .7  A,  depending  on  the  basis,  there  is 
one  imaginary  frequency.  Displacement  along  this  mode  corresponds  to  an  asymmetric  axial 
stretch  leading  to  dissociation  into  the  fluorine  anion  (F")  and  the  oxypentafluoride  cation  (OF5^) 
in  C4V  symmetry.  This  conformation  for  OF5'‘'  is  a  saddle  corresponding  to  the  Berry 
pseudorotation  state  in  the  D311  C4V  D3J1  ligand  interchange  process.  The  calculated  energy 
barrier  for  the  pseudorotation  process  in  OF5''’  is  6.5  kJ/mol  at  the  HF/6-3 1  G(d)  level  of  theory. 
This  is  similar  to  the  calculated  barrier  of  8.8  kJ/mol  for  the  corresponding  process  in  NF5.  At 
the  HF  level  of  theory,  the  covalent  OF6  structure  is  energetically  favored.  The  predicted  long 
bond  length  at  the  MP2  level  of  theory,  using  diffuse  basis  flmctions,  results  in  a  reference  state 
which  is  UHF  unstable.  Thus,  MP2  may  not  provide  an  adequate  description  of  the  electronic 
structure  of  this  system  where  there  appears  to  be  a  stability  competition  between  the  covalent 
D3h  structure  and  the  0F5''T‘  ion-pair.  The  stability  of  OFg  is  therefore  similar  to  that  for  NF5. 

A  gas  phase  species  may  exist,  but  the  ionic  form  appears  to  be  favored  in  the  solid. 
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Table  1.  Gradient  Optimized  Structures  and  Energies  for  OF5 


Theory 

R(A) 

HF/6-31G 

1.6126 

HF/6-31G(d) 

1.5407 

HF/6-31G(2df) 

1.5293 

HF/6-31+G(d) 

1.5457 

HF/6-311G(d) 

1.5470 

MP2/6-31G(d) 

1.6673 

MP2/6-31+G(d) 

1.6956 

MP2/6-311G(d) 

1.7210 

MP2/6-3ll+G(d) 

1.7452 

Efh')  ZPEOi) 


-670.2323  325 

0.0168  625 

-670.3343  435 

0.0183  033 

-670.3586  884 

-670.3605  824 

0.0177  040 

-670.5047  114 

0.0166  128 

-671.8802  131 

0.0128  721 

-671.9594  044 

— 

-672.4353  062 

* 

-672.4817  277 

*  1  imaginary  frequency  -  asymmetric  axial  stretch;  —  not  calculated 


3.2  Azide-like  Structures 

3.2.1  N(N3)3,  IIN(N3)2,  N(N3)2'  and  N(N3)/.  The  successful  synthesis  of  halogen 
azides  and  dinitramide  structures  has  led  to  an  examination  of  the  possibility  of  stability  of  other 
unusual  nitrogen  structures.  Previous  theoretical  studies  have  examined  the  stability  of  several 
energetic  forms  of  Nj^.  Ng  has  been  shown  to  be  stable  as  a  hexazabenzene  structure  at  the  SCF 
level  of  theory  (Refs.  1-3)  but  subsequent  calculations  at  a  correlated  level  of  theory  indicate  this 
conformation  to  be  a  second-order  saddle  point  (Refs.  1, 4).  Several  studies  of  Ng  as  an  open 
chain  structure  have  been  reported  (Refs.  5, 6).  The  lowest  conformation  appears  to  have  C2 
symmetry  and  exists  as  a  weak  dimer  of  the  azide  radical.  Engelke  (Ref.  2)  has  examined  several 
other  conformations  for  Ng;  all  appear  to  lie  higher  in  energy  than  the  open  chain  diazide  form. 
An  experimental  study  by  Vogler,  et  al.  (Ref.  7)  postulated  that  Ng  is  formed  as  a  product  of  an 
elimination  reaction  of  a  diazido  platinum  complex  at  low  temperatures.  No  direct  evidence  for 
formation  of  a  stable  form  of  Ng  has  been  reported.  Studies  of  other  polyhedral  forms  of  Nx 
have  been  reported  (Refs.  4, 8-12).  Of  particular  interest  is  the  tetrahedral  N4  structure,  since  the 
corresponding  P4  structure  is  more  stable  than  P2.  Recent  calculations  (Ref.  4)  at  a  correlated 
(MBPT-2)  level  of  theory  suggest  a  barrier  of  approximately  50  kcal/mol  for  decomposition  of 
N4.  Octaazacubane  [Ng]  has  also  been  examined  at  both  the  SCF  and  correlated  levels  of  theory 
(Refs.  4,  8-10).  This  highly  strained  molecule  is  predicted  to  have  a  gas  phase  heat  of  formation 
of  -  470  kcal/mol  (Ref  4)  and  decomposes  to  4  N2  via  a  high  activation  barrier.  Unlike  cubane, 
however,  Ng  (Oh  symmetry)  has  no  convenient  intermediate  stmctores  such  as  those  found  by 
Eaton,  et  al.  (Refs.  1 1, 12)  for  cubane  synthesis. 

In  contrast  to  the  numerous  studies  of  even-membered  nitrogen  compounds,  little  has  been 
reported  on  odd-membered  nitrogen  structures  beyond  the  azide  radical  [N3].  Nguyen,  et  al. 
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(Refs.  13, 14)  have  examined  N5'  as  an  open  chain  fomi,  and  as  pentazene  and  pentazole 
stractures.  The  lowest  stable  conformation  appears  to  be  the  open  chain  form,  which  is  a  weak 
van  der  Waals  complex  of  N2  +  N3'.  The  pentazole  anion  appears  stable  at  the  RHF  level  of 
theory,  lying  some  16  kcal/mol  above  the  van  der  Waals  structure.  Further  studies  at  a  corrected 
level  of  theory  are  needed  to  better  define  the  relative  energetics  of  these  N5‘  conformations. 

The  next  likely  candidate  for  stability  among  odd-membered  nitrogen  compounds  is  the 
diazidamine  anion  [N7'].  This  anion  could  exist  as  a  N(N-N2)2’  structure,  similar  to  the 
dinitramide  anion  [N(N02)2’]>  or  as  a  N(N3)2”  structure  with  C2  symmetry  formed  by  the  linkage 
of  two  azide  groups  to  a  central  nitrogen  atom.  Therefore  ab  initio  calculations,  including 
electron  correlation,  were  carried  out,  for  both  possible  conformations.  As  shown  in  the  next 
paragraph,  the  C2  structure  proved  to  be  a  stable  minimum.  UTRC  then  examined  protonation  of 
this  structure  to  form  hydrogen  diazidamine  [HN(N3)2],  the  simplest  stable  molecule  of  this 
anion.  Finally,  the  triazidamine  molecule  [N(N3)3]  was  examined  and  a  stable  C3  structure  was 
formd  for  this  highly  energetic  molecule. 

Ab  initio  calculations  were  carried  out  for  Ny"  in  both  C2  and  Cg  symmetries  using  the 
Gaussian  92  programs  (Ref  15).  Full  geometry  optimization  at  the  RHF/6-31G(d)  level  of 
theory,  followed  by  a  calculation  of  harmonic  firequencies,  indicated  that  the  N(N-N2)2'  structure 
was  not  a  vibrationally  stable  structure  (two  imaginary  vibrational  firequencies).  However,  at  this 
same  level  of  theory,  the  open  chain  diazidamine  form  of  Ny’  was  found  to  be  stable  in  C2 
symmetry,  exhibiting  all  real  harmonic  firequencies.  Additional  calculations,  improving  the  basis 
to  include  diffuse  fimctions  and  geometry  optimization  at  the  MP2  level  of  theory,  were  then 
carried  out.  The  results  are  given  in  Tables  2  and  3,  respectively,  and  the  stable  C2  conformation 
is  illustrated  in  Figure  1 .  At  the  MP2  correlated  level  of  theory,  the  central  N1-N2  bond  distance 
(1.45  A)  is  slightly  lengthened  (0.08  A)  at  MP2[6-31+G(d)]  over  that  previously  reported  for  the 
dinitramide  anion.  Other  known  compormds  of  nitrogen  exhibit  even  larger  N-N  bond  lengths: 
N2O3  (1.86  A),  N2O2  (1.78  A),  and  N2H4  (1.47  A).  The  N2-N3  and  N3-N4  bond  lengths  are 
typical  of  those  found  in  hydrogen  azide  and  several  halogen  azide  compormds  (Refs.  16, 17). 

The  central  bond  angle  is  close  to  that  for  sp^  hybridization  on  the  central  nitrogen  atom.  A 
CHELPG  charge  analysis  (Ref  18)  indicates  an  excess  (0.41)  negative  charge  on  the  central 
nitrogen  atom. 

Previous  studies  have  indicated  a  rich  structural  chemistry  for  hydrogen  dinitramide 
(Ref  19).  Accordingly,  UTRC  examined  the  stability  of  the  simplest  diazidamine  compound, 
hydrogen  diazidamine.  Geometry  optimization  was  carried  out  at  the  RHF/6-3  lG(d)  and 
MP2/6-31G(d)  levels  of  theory.  The  vibrational  analysis  reveals  the  stable  Cg  structure  given  in 
Table  4  and  illustrated  in  Figure  2.  We  find  geometric  parameters  for  HN(N3)2  that  are  very 
similar  to  those  shown  in  Table  2  for  the  diazidamide  anion  [N(N3)2']. 

In  a  parallel  fashion,  UTRC  examined  the  structural  stability  of  triazidamine,  a  compound 
analogous  to  trinitramine  [N(N02)3]  that  was  recently  predicted  structurally  to  be  at  least  as 
stable  as  dinitrogen  pentoxide  [N2O5  ].  Ab  initio  geometry  optimization  calculations  of  N(N3)3 
were  first  carried  out  at  the  RHF/6-3  lG(d)  level  of  theory.  Harmonic  frequency  calculations 
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Table  2.  Calculated  Enei^es  and  Geometries  of  the  Diazidamide  Anion  IN(N3)2’] 


Level  of  Theory 

RHF/6-31Gfd) 

RHF/6-31+Gfdl 

MP2/6-31G('d) 

MP2/6-31+G('d') 

Energy 

-380.866  526 

-380.893  845 

-382.049  745 

-382.097  251 

Geometry 

Ri2 

1.4180 

1.4089 

1.4529 

1.4425 

^^23 

1.1974 

1.1976 

1.2465 

1.2472 

R34 

1.1329 

1.1310 

1.1754 

1.1743 

®2'12 

101.81 

102.70 

100.58 

101.21 

®123 

113.18 

112.82 

115.69 

113.65 

®234 

176.56 

177.08 

170.02 

171.50 

dihri?^ 

-180.00 

-173.72 

-150.11 

-156.43 

dihi234  -180.00  -154.69  -143.95 

Energies  are  in  hartrees  and  geometries  are  in  angstroms  and  degrees. 

RHF/6-31G(d)  zero-point  energy  (scaled)  =  0.026488. 

-137.65 

Table  3.  Calculated  Harmonic  Frequencies  for  the  Diazidamide  Anion  [N(N3)2"] 

Level  of  Theory 

RHF/6-31G('d') 

RHF/6-31+G('d') 

MP2/6-31G('d) 

MP2/6-31+G('d') 

Freauencies  FIR  intensitvl 

(vibrational  frequencies  are  in  cm“^  and  IR  intensities  are  in  km/mol) 

ffll  (B) 

56  [5] 

74  [8] 

35  [3] 

44  [6] 

®2  (A) 

72  [0] 

107  [0] 

63  [0] 

85  [0] 

<»3  (A) 

140  [1] 

145  [0] 

150  [1] 

145  [0] 

(04  (B) 

278  [24] 

282  [23] 

218  [25] 

227  [23] 

<05  (A) 

459  [4] 

465  [5] 

396  [0] 

379  [0] 

<06  (B) 

526  [33] 

492  [30] 

385  [12] 

326  [13] 

<07  (A) 

523  [0] 

510  [0] 

425  [0] 

408  [0] 

<»8  (A) 

622  [4] 

629  [2] 

610  [7] 

593  [2] 

©9  (B) 

799  [118] 

796  [159] 

743  [60] 

737  [4] 

©10  (A) 

991  [20] 

1001  [16] 

914  [4] 

911  [1] 

©11  (B) 

982  [114] 

1024  [173] 

869  [47] 

899  [78] 

©12  (B) 

1421  [7] 

1428  [9] 

1162  [2] 

1159  [2] 

©13  (A) 

1473  [7] 

1476  [14] 

1188  [9] 

1184  [21] 

©14  (B) 

2318  [2634] 

2304  [2970] 

2280  [1197] 

2255  [1422] 

©15  (A) 

2359  [190] 

2388  [256] 

2291  [14] 

2263  [47] 
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Figure  1 

Diazidamide  Anion  [N(N3)2~1  -  Cj  Symmetry  at  the  MP2/6-31G(d)  Level  of  Theory 


Previous  studies  have  indicated  a  rich  structural  chemistry  for  hydrogen  dinitramide 
(Ref  19).  Accordingly,  UTRC  examined  the  stability  of  the  simplest  diazidamine  compound, 
hydrogen  diazidamine.  Geometry  optimization  was  carried  out  at  the  RHF/6-3  lG(d)  and 
MP2/6-31G(d)  levels  of  theory.  The  vibrational  analysis  reveals  the  stable  Cg  structure  given  in 
Table  4  and  illustrated  in  Figure  2.  We  find  geometric  parameters  for  HN(N3)2  that  are  very 
similar  to  those  shown  in  Table  2  for  the  diazidamide  anion  [N(N3)2^. 


Table  4.  Calculated  Energies  and  Geometries  of  the  Diazidamine  Molecule  [HN(N3)2] 


Level  of  Theory 

RHF/6-3  IGtd) 

MP2/6-31Gtd) 

Energy  (hartrees) 

-381.471  597 

-382.633  146 

Geometry  (angstroms  and  deerees^l 

Ri2 

1.4249 

1.4665 

Ri5 

0.9989 

1.0213 

%3 

1.2508 

1.2679 

1.0971 

1.1563 

®215 

104.25 

101.16 

®123 

108.14 

108.39 

®234 

173.84 

171.29 

d‘1'2’152 

114.89 

111.10 

dih5j23 

-136.24 

-153.03 

dihj234 

-178.19 

-173.05 

RHF/6-3  lG(d)  zero-point  energy  (scaled)  =  0.041008. 
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Figure  2 

Diazidamine  Molecule  [HN(N3)2l  -  Cj  Symmetry  at  the  MP2/6-31G(d)  Level  of  Theory 


In  a  parallel  fashion,  UTRC  examined  the  structural  stability  of  triazidamine,  a  compound 
analogous  to  trinitramine  [N(N02)3]  that  was  recently  predicted  structurally  to  be  at  least  as 
stable  as  dinitrogen  pentoxide  [N2^5  initio  geometry  optimization  calculations  of  N(N3)3 
were  first  carried  out  at  the  RHF/6-3  lG(d)  level  of  theoiy.  Harmonic  frequency  calculations 
confirmed  a  vibrationally  stable  C3  structure,  as  illustrated  in  Figure  3. 


Triazidamine  Molecule  [N(N3)3|  -  C3  Symmetry  at  the  RHF/6-31G(d)  Level  of  Theory 
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An  improved  geometry  was  optimized  at  the  MP2/6-3  lG(d)  level,  and  again  the  C3  structure  of 
N(N3)3  was  found  to  be  vibrationally  stable.  Inclusion  of  diffuse  functions  at  the  RHF  level 
resulted  in  only  minor  changes  in  predicted  bond  lengths  and  angles.  The  results  are  summarized 
in  Tables  5  and  6. 

Table  5.  Calculated  Energies  and  Geometries  of  the  Triazidamine  Molecule  [N(N3)3l 


Level  of  Theory 

RHF/6-3  IGld! 

RHF/6-3 1+Gfd! 

MP2/6-31Gfd! 

Energy  rtjartreesl 

-544.061  595 

-544.128  121 

-545.773  964 

GeometryCangstroms  and  degrees! 

Ri2 

1.4219 

1.4206 

1.4626 

%3 

1.2514 

1.2515 

1.2713 

R34 

1.0966 

1.0961 

1.1551 

®2'12 

106.57 

106.73 

103.87 

®123 

107.41 

107.46 

106.46 

®234 

174.02 

173.94 

171.55 

dih2<i23 

147.61 

150.49 

147.30 

dihi234 

174.65 

174.19 

177.34 

RHF/6-3  lG(d)  zero-point  energy  (scaled)  =  0.042438. 


As  in  ammonia  (NH3),  nitrogen  trifluoride  (NF3),  trimethylamine  and  trinitramine 
(Ref  20),  UTRC's  calculations  show  that  triazidamine  is  pyramidal,  not  planar.  The  N-N  central 
bond  distance  is  1.46  A  in  N(N3)3,  very  close  to  that  found  in  hydrazine  (1.47  A)  and  in  the 
diazidamine  anion  (1.44  A),  and  considerably  shorter  than  the  corresponding  N-N  bond  distance 
in  trinitramine  (1 .54  A).  This  suggests  weak  electrostatic  interaction  between  the  azide  groups  in 
N(N3)3.  In  covalent  azides,  such  as  hydrogen  azide  (HN3)  and  fluorine  azide  (FN3),  thermal 
decomposition  occurs  by  cleavage  of  the  weak  N-N  double  bond: 


HN3  N2  +  NH 

(1) 

FN3  —>  N2  +  NF 

(2) 

The  barrier  to  decomposition  is  small  for  HN3  and  ~10.9  kcal/mol  for  FN3,  much  less  than  the  N- 
H  and  N-F  bond  strengths.  In  HN(N3)2  and  N(N3)3,  the  single  bond  to  the  central  nitrogen  is 
probably  much  weaker  (-40  kcal/mol  in  nitramides),  but  still  larger  than  the  barrier  for: 


HN(N3)2  ^ 

N2  +  fragments 

(3) 

N(N3)3  ^ 

N2  +  fragments 

(4) 
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Table  6.  Calculated  Harmonic  Frequencies  for  the  Triazidamine  Molecule  [N(N3)3] 

Level  of  Theory  RHF/6-31G(d)  RHF/6-31G(d)  MP2/6-31GIdl 


and  IR  intensities  are  in  km/mol) 


©i  (E) 
©2  (E) 

©3  (A) 

©4  (A) 

©5  (E) 

cog  (E) 

©7  (E) 

©g  (E) 

©9  (A) 

©10  (A) 
©11  (E) 

©12  (E) 

©13  (A) 

©14  (E) 

©15  (E) 

®16  (E) 

©17  (E) 

©Ig  (A) 

©19  (E) 

©20  (E) 
®21 

©22  (E) 
©23  (E) 

©24  (A) 


36  [1]  39 

36  [1]  39 

52  [0]  59 

155  10]  164 

252  [2]  251 

252  [2]  251 

468  [12]  466 

468  [12]  466 

536  [15]  539 

593  [0]  592 

644  [16]  626 

644  [16]  626 

649  [1]  635 

796  [13]  793 

796  [13]  793 

1130  [23]  1129 

1130  [23]  1129 

1109  [4]  1098 

1217  [314]  1220 

1217  [314]  1220 

1276  [107]  1272 

2457  [739]  2454 

2457  [739]  2454 

2490  [403]  2485 


11] 

19 

[0] 

in 

19 

[0] 

[1] 

35 

[0] 

[0] 

129 

[0] 

[2] 

212 

[1] 

[2] 

212 

[1] 

[11] 

413 

[6] 

[11] 

413 

[6] 

[15] 

486 

[3] 

10] 

524 

[0] 

[13] 

526 

[2] 

[13] 

526 

[2] 

[0] 

549 

[0] 

[12] 

713 

[24] 

[12] 

713 

[24] 

[39] 

897 

[83] 

[39] 

897 

[83] 

[5] 

931 

[0] 

[321] 

1171 

[94] 

[321] 

1171 

[94] 

[126] 

1221 

[15] 

[771] 

2380 

[375] 

[771] 

2380 

[375] 

[445] 

2387 

[103] 

Previous  calculations  on  FN3  have  shown  that  the  experimentally  observed  barrier  (11.5  kcal/mol) 
for  reaction  (2)  is  well  reproduced  at  the  RHF/6-3  lG(d)  level  of  theory.  Gradient  optimized 
transition  states  for  reactions  (3)  and  (4)  have  been  located  at  the  RHF/6-3  lG(d)  level,  and  the 
barriers  to  the  decomposition  processes  (3)  and  (4)  are  9.9  and  10.8  kcal/mol,  respectively.  The 
preliminary  conclusion  is  that  the  shock  sensitivity  of  HN(N3)2  and  N(N3)3  may  be  comparable  to 
that  of  fluorine  azide. 

To  complete  the  structural  possibilities  of  polyazamides,  the  tetraazidammonium  cation 
[N(N3)4‘^]  was  also  examined.  Ab  initio  geometry  calculations  were  carried  out  at  the  RHF/6- 
3  lG(d)  level  of  theory.  A  harmonic  frequency  analysis  indicated  a  vibrationally  stable  D2d 
structure,  as  illustrated  in  Figure  4.  An  improved  optimized  geometry  was  then  found  at  the 
MP2/6-31G(d)  level  of  theory.  The  calculated  geometries  and  harmonic  frequencies  are  shown 
in  Tables  7  and  8,  respectively. 
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Tetraazidammonium  Cation  (N(N3)4''’]  -  D2(i  Symmetry  at  MP2/6-31G(d)  Level  of  Theory 


Table  7.  Calculated  Energies  and  Geometries  of  the 
Tetraazidammonium  Cation  [N(N3)4^] 


Level  of  Theory 

RHF/6-31G(d) 

MP2/6-31G(d) 

Energy  Thartrees) 

-707,070  221 

-709.216  858 

Geometry  tangstroms  and  degrees) 

Ri2 

1.4359 

1.4651 

%3 

1.2929 

1.2919 

R34 

1.0850 

1.1489 

®2'12 

100.08 

99.15 

0123 

106.81 

106.05 

®234 

171.75 

169.88 

RHF/6-31G(d)  zero-point  energy  (scaled)  =  0.055930. 


Table  8.  Calculated  Harmonic  Frequencies  for  the  Tetraazidammonium  Cation  [N(N3)4'''] 


RHF/6-31Gfdt  Level  of  Theory 

Al:  113  10),  348  [OJ,  547  |0].  943  (O),  1112  (0),  2599  [0] 

A2;  118  [0],  582  |0| 

Bl;  82  10),  423  |0],  614  |0] 

B2:  163  |0].  443  16],  592  11],  1074  126],  1231  11],  2589  152] 

E:  87  12],  226  13],  473  111],  615  139],  842  149],  965  196],  1161  ]  231],  2589  1292] 


Vibrational  frequencies  are  in  cm*'  and  IR  intensities  ]  ]  are  in  km/mol. 
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Comparing  the  geometries  of  the  tetraazidammonium  cation  with  that  found  for  the  neutral 
triazidamine  molecule,  we  find  nearly  identical  N-N  bond  lengths  to  the  central  nitrogen  atom  at 
both  the  RHF  and  MP2  level  of  theory.  This  indicates  that  there  is  little  additional  electrostatic 
repulsion  (crowding)  resulting  from  the  added  N3''^  group  in  the  cation.  Similarly,  the  central  and 
end  nitrogen  atom  bond  lengths  in  the  azide  groups  do  not  significantly  differ  between  the  neutral 
triazidamine  and  tetraazidammonium  cation  stmcture. 

The  heat  of  formation  of  HN(N3)2  can  be  extracted  fi’om  the  enthalpy  change  for  the 
isogyric  reaction: 

NH3  +  2HN3  ^  HN(N3)2  +  2H2  (5) 

Using  a  calculated  heat  of  reaction  of  64.6  kcal/mol  at  the  MP2/6-3  lG{d)  level  of  theory  for 
reaction  (5)  and  the  experimental  heats  of  formation  of  NH3  and  HN3  (Ref  21),  AHf°  (0  K)  = 
+198.9  kcal/mol  for  HN(N3)3. 

The  heat  of  formation  of  N(N3)3  can  be  estimated  from  the  heat  of  reaction  for 

5  N2  Nio  (6) 

or,  alternatively,  from  the  isogyric  reaction 

NH3  +  3  HN3  ^  N(N3)3  +  3  H2  (7) 

For  reaction  (7),  the  heat  of  formation  of  triazidamine  was  estimated  from  the  calculated  MP2/6- 
3  lG(d)  heat  of  reaction  and  the  experimental  heats  of  formation  of  HN3  and  NH3  (Ref  21).  The 
results  are  summarized  in  Table  9,  which  indicates  that  N(N3)3  represents  a  very  energetic  and 
possibly  explosive  material,  AHf®  (N(N3)3)  =  +298  kcal/mol.  This  energy  content  arises, 
however,  fi’om  the  large  difference  in  the  N-N  single  and  double  strengths  compared  to  that  of  the 
triple  bond,  rather  than  as  the  strain  energy  that  arises  in  a  structure  such  as  octaazacubane. 
Further,  the  building  block  for  this  molecule  is  the  diazidamine  anion,  which  may  possibly  be 
synthesized  starting  from  dimethylamine. 


HN(N3)2 

N(N3)3 


Table  9.  Calculated  Heat  of  Formation  for  the 


Diazidamine  [HN(N3)2]  and  Triazidamine  [N(N3)3]  Molecules 


Reaction 


MP2/6-31G(d)  Level  of  Theoiy 

Equation 


AHf°fOK) 


NH3  +  2HN3  ^  HN(N3)2  +  2H2 

NH3  +  3  HN3  N(N3)3  +  3  H2 


(5)  198.9  kcal/mol 

(7)  298.1  kcal/mol 
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3.2.2  NONj  and  NO2N3.  Calculations  on  several  nitro  azides  were  carried  out.  Klapotke 
and  coworkers  (Ref.  22)  have  recently  reported  the  synthesis  of  nitrosyl  azide  [NO^N3’].  This 
energetic  compound  is  stable  at  low  temperatures,  existing  as  a  chain-type  structure  rather  than  as 
a  cyclic  form.  Calculations  for  NO^N3’  and  for  N02^N3'  (nitryl  azide)  were  carried  out  at  the 
MP2/6-3 1+G(d)  level  of  theory.  The  optimized  structures  (both  planar)  are  shown  in  Figures  5 
and  6,  respectively. 


It  is  of  interest  to  compare  the  relative  value  of  several  oxidizers,  including  the  nitro  azides, 
in  improving  the  performance  characteristics  of  solid  propellant  formulations.  Heats  of  formation 
for  all  of  these  compounds  have  been  calculated  at  the  MP2/6-3 1+G(d)  level  of  theory  from 
isodesmic  reactions.  The  use  of  isodesmic  reactions  permits  cancellation  of  errors  that  can 
otherwise  accumulate  if  heats  of  formation  are  based  on  atomization  energies.  Solid  densities 
were  estimated  using  Cady's  method  (Ref  23)  for  evaluating  explosives.  These  calculated  data 
are  shown  in  Table  10. 
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Table  10.  Thermochemistry  of  Advanced  Propellants 


Compound 

Triprismane 

Nitrotriprismane 

Dinitrotriprisniane 

Tetranitrohydrazine 

Trinitiamine 

Diazidamine 

Triazidamine 

Diazidoborane 

Imidazidoborane 

Nitrosyl  azide 

Nitiyl  azide 

3-nitro-l,2,4-triazol-5-one  (NTO) 
l-anuno-3-nitro-l,2,4-triazol 


Formula 

Mol.  Wt. 

W 

78.12 

CgHjNOj 

123.11 

C6H4(N02)2 

168.11 

N2(N02)4 

212.04 

N(N02)3 

152.02 

HN(N3)2 

99.06 

N(N3)3 

140.07 

HB(N3)2 

95.86 

HNBN3 

67.85 

NON3 

72.03 

NO2N3 

88.03 

C2H2N4O3 

130.06 

C2H3N5O2 

129.08 

Density 

em/cm^ 

AH£(25<»C) 

kcal/mol 

1.29 

+134.7 

1.55 

+125.3 

1.71 

+121.5 

2.17 

+140.0 

2.15 

+90.0 

1.56 

+199.0 

1.58 

+302.0 

1.52 

+108.0 

1.34 

+56.0 

1.54 

+105.0 

1.76 

+90.0 

1.75 

-3.2 

1.61 

+55.6 

Density  is  calculated  by  the  method  of  Cady  (Ref.  23). 


The  new  energetic  materials  are  separated  into  three  groups:  oxidizers,  monopropellants, 
and  oxygen-deficient  energetic  materials.  These  compounds  are  formulated  in  propellants  and  the 
Igp's  of  these  propellants  are  calculated  through  the  AFPL  theoretical  Igp  program.  The  oxidizers 
are  calculated  as  ammonium  dinitramide  (ADN)  substitutes  in  an  ammonium  dinitramide/glycidyl 
azide  polymer  (ADN/GAP)  smokeless  propellant.  The  monopropellants  replace  HMX  in  a 
smokeless/PEG  (polyethylene  glycol)  propellant.  The  oxygen-deficient  materials  are  formulated 
in  an  aluminized/HTPB  (hydroxyl  terminated  polybutadiene)  propellant  as  replacements  for  HMX. 

Oxidizers.  Tetranitrohydrazine  [N2(N02)4],  trinitramine,  nitrosyl  azide,  and  nitryl  azide  are 
the  oxidizers  formulated  in  a  GAP/smokeless  propellant.  The  calculated  Igp's  are  shown  in 
Table  1 1.  An  average  of  30  to  40  sec  improvement  in  the  specific  impulse  over  the  formulation 
containing  ADN  as  an  oxidizer  is  achieved.  Similar  improvement  is  also  seen  in  density  Igp.  The 
large  improvement  in  performance  also  comes  with  much  lower  percentages  of  oxidizer  contents. 
Lower  solid  contents  make  propellant  processing  easier.  Lower  solid  contents  also  produce 
propellants  with  better  physical  properties. 

Monopropellants.  Triazidamine  (TAA),  a  monopropellant,  is  used  in  replacement  of 
HMX,  also  a  monopropellant,  in  a  selected  formulation.  A  generic  energetic  smokeless/PEG 
formulation  is  chosen  for  comparison.  The  base  formula  is  not  optimized.  This  base  formulation 
contains  HMX  (60  %),  PEG  (13  %),  NG  (nitroglycerin)  (25  %),  and  N-3200  (2  %).  As  shown  in 
Table  1 1,  a  50  sec  improvement  in  Igp  is  achieved. 
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Oxygen-Deficient  Energetic  Materials.  Diazidoborane  [HB(N3)2],  imidazidoborane 
[HNBN3],  l-ainino-3-nitro-l,2,4-triazol  [C2H3N5O2],  and  the  two  nitrotriprismanes  are  oxygen- 
deficient  energetic  materials.  A  good  way  has  not  yet  been  found  to  incorporate  these  materials 
into  existing  propellant  formulations  for  a  quick  performance  comparison.  For  now,  an 
aluminized/HTPB  system  is  chosen  as  a  base  propellant.  The  propellant  has  the  following 
formulation:  HTPB  (8.25  %),  DDI  (dimeryl  diisocyanate)  (1.25  %),  aluminum  (20  %),  AP 
(58  %),  and  HMX  (12  %).  Again,  this  generic  formulation  is  used  for  comparison  only  and  is  not 
optimized.  Diaadoborane  and  l-amino-3-nitro-l,2,4-triazol  are  used  to  replace  HMX  in  the  base 
formula.  The  percentages  of  aluminum  and  AP  are  adjusted  to  maximize  the  calculated 
performance.  Table  1 1  lists  the  calculated  results.  UTRC  believes  that  an  improvement  in 
performance  can  be  realized  upon  finding  a  suitable  propellant. 


Table  11.  Theoretical  Performance  Characteristics  (Ijp)  of  Advanced  Propellants 

Oxidizers 


Ingredient 

Isp,  sec 

Isp  (vac) 

Isp*Density  Isp(vac)*Density 

Density 

%Oxidi: 

AP 

255 

277 

16.2 

17.5 

1.75 

78 

ADN 

265 

287 

16.2 

17.6 

1.69 

83 

N2(N02)4 

293 

320 

18.6 

20.3 

1.75 

62 

N(N02)3 

293 

319 

18.2 

19.9 

1.71 

64 

NON3 

305 

331 

16.4 

17.8 

1.48 

81 

NO2N3 

294 

319 

17.1 

18.6 

1.61 

74 

MononroDellants 

Ingredient 

Isp,sec 

Isp  (vac) 

Isp*Density 

Isp  (vac)*Density 

Density 

HMX  (base) 

243 

261 

14.7 

15.8 

1.67 

TAA 

296 

318 

16.2 

17.4 

1.51 

Oxveen-Deficient  Energetic  Materials 

Ingredient 

Isp,  sec 

Isp  (vac) 

Isp*Densily 

Isp  (vac)*Density 

Density 

HMX  (base) 

266 

291 

17.7 

19.3 

1.84 

ANT 

265 

291 

17.2 

18.9 

1.79 

DAB 

266 

291 

17.0 

18.7 

1.76 
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3  23  ON(N,)j^  Calculations  were  utoted  on  the  oxodiaidanude^^^ 

fo.«suggT^nhyIeffBottaro(SRI).hat^^^ 

vibrationaUy  stable  smtchne  has  been  locate  at  boft  ^OT/6  M 

of  theory.  This  structure  exhibits  C2  symmeiy  si  vibrational  fiequendes 

^  14F/6-'?  1  level  is  shown  in  Table  12.  i  ne 

Table  13.  A  pictorial  representation  of  this  new  cation  is  shown  inFigure 

Table  12.  Calculated  Geomrtry  uf  the  Oxodiazidamide  Cation  lONfNsh^ 


Oi 


-N2; 


■Ns' 

•N3. 


.N4 


N4. 


N5 


N5. 


Geometiy  is  in  angstroms  and  degrees 


Ri2 

^23 

R34 

R45 

®123 

®234 

6345 

1.150 

1.365 

1.305 

1.083 

117.54 

110.37 

169.17 

X1234  '^2345 

142.51  188.58 


Table  13.  Calculated  Harmonic  Frequmtcics  mtd  Infrared  luteusities  for  ON(Nsh" 


Frequencies  are  in  cm 


fi  and  IR  intensities  [  ]  are  in  km/mol 


Mode  Freguencx 

B  71.2  [2] 

230.7  13] 

545.0  [8] 

572.2  [8] 
626.5  [1] 
874.1  [46] 

1042.8  [53] 

1285.7  [264] 

2626.9  [74] 


Mode  Frequencx 

A  87.7  [2] 

189.8  [3] 

459.7  [0] 

561.9  [30] 

745.8  [1] 
881.6  [1] 

1118.1  [264] 

2005.5  [380] 
2636.0  [94] 
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Figure  7 

Optimized  Structure  of  the  Oxodiazidamide  Cation 


|on(N3)2+i 


3.2.4  Azidoboranes.  Ab  initio  calculations  were  carried  out  on  several  azidoborane 
structures.  Wiberg  and  Michaud  (Ref.  24)  have  shown  that  HN3  and  diborane  (B2H5)  react  in 
ether  by  abstracting  H2  to  form  azidoboranes.  The  gas  phase  reaction  could  proceed  as 

H 

N3  /  \  N3 


\  / 

\  / 

B2H6  +  4HN3  ^  B 

—  B 

+ 

2H2 

(8) 

/  \ 

/  \ 

N3  \  /  N3 

H 

provided  the  reaction  products  are  stable.  Ab  initio  calculations  were  performed  to  explore  the 
possible  existence  and  stability  of  the  tetraazidoborane  (TADB)  structure  suggested  in  reaction 
(8).  The  RHF/6-3  lG(d)  optimized  geometry  for  TADB  is  shown  in  Table  14,  and  the  calculated 
vibrational  frequencies  and  intensities,  which  can  be  used  to  characterize  experimental  absorption 
bands,  are  given  in  Table  15.  The  calculations  indicate  the  likely  existence  of  a  vibrationally  stable 
D2  structure  for  TADB.  Calculated  heats  of  reaction  and  heats  of  formation  are  reported  in 
Table  16.  These  structures  include  vibrational  zero-point  corrections  obtained  from  the  calculated 
vibrational  frequencies  (scaled  by  0.8929).  The  MP2/6-3  lG(d)//RHF/6-3  lG(d)  heat  of  formation 
of  TADB,  given  in  Table  16,  indicates  that  reaction  (8)  is  exothermic  by  73.2  kcal/mol. 
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Table  14.  Calculated  Geometries  for  Several  Azidoborane  Structures 


RHF/6-31G(d)  Level  of  Theoiy 
Geometries  are  in  angstroms  and  degrees 


Tetraazidoborane 


TADB  (P2) 

B-B 

1.821 

H-B 

1.332 

B-Nj 

1.464 

N1-N2 

1.230 

N2-N3 

1.097 

H-B-Ni 

103.5 

B-N1-N2 

118.8 

NPN2-N3 

174.7 

H-B-N1-N2 

-168.2 

B-N1-N2-N3 

-179.5 

Diazidoborane 


DAB  (C2v) 

H-B 

1.179 

B-Ni 

1.439 

N1-N2 

1.237 

N2-N3 

1.095 

H-B-Nj 

115.4 

B-N1-N2 

123.6 

N1-N2-N3 

173.1 

Imidazidoborane 
HNBN3  (Cg) 


H-N 

0.980 

N-B 

1.225 

B-Ni 

1.409 

N1-N2 

1.233 

N2-N3 

1.094 

H-N-B 

180.0 

N-B-H| 

177.4 

B-N1-N2 

120.1 

N1-N2-N3 

173.9 

N-B-N1-N2 

180.0 

Table  15.  Calculated  Vibrational  Frequencies  and 
Intensities  for  Several  Azidoborane  Structures 

RHF/6-31(d)  Level  of  Theoiy 
Frequencies  are  in  cm'*  and  intensities  [  ]  are  in  km/mol 


DAB  Diazidoborane 

Ai  28451162],  2588[924],  1401[556],  919[16],  824[81],  506[10],  122(0] 

A2  674(0],  60(0] 

Bj  952(59],  686(63],  311(3] 

B2  2532(194],  1435(900],  1279(13],  1043(27],  693(1],  181(0] 

TADB  Tetraazidohnrane 

A  2580(0],  2260(0],  1442(0],  1103(0],  835(0],  674(0],  532(0],  2622(0],  186(0],  79(0],  42(0] 

Bi  2533(633],  2003(106],  1442(262],  1158(18],  730(47],  668(33],  497(25],  315(1],  155(3],  34(1] 

B2  2537(440],  1744(4],  1450(843],  1312(385],  1080(90],  724(1],  667(8],  453(0],  181(0],  110(1],  33(0] 

B3  2560(2145],  1831(2335],  1357(3117],  1160(391],  877(659],  800(368],  668(90],  485(15],  134(10],  48(2] 
HNBN3  Imidazidoborane 

A'  4146(336],  2592(728],  2178(882],  1432(300],  900(30],  695(11],  559(154],  517(19],  153(0] 

A"  651(19],  490(0],  435(210] 
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Table  16.  Calculated  Energies  and  Thermochemistry  for  Several  Azidoborane  Structures 


Level  of  Theory 

E(hartrees) 

AHf  (kcal/mol) 

AHf  (0  K)(kcal/mol) 

TABD  Tetraazidoborane  tSee  Ea.  8) 

RHF/6-31G(d) 

-703.68903 

-33.3 

266.2 

MP2/6-31G(d)//RHF/6-31G(d) 

-705.83842 

-75.0 

224.5 

MP2/6-3  l+G(d)//RHF/6-3  lG(d) 

-705.87845 

-73.2 

226.3 

DAB  Diazidoborane  fSee  Ea.  91 

RHF/6-3  lG(d) 

-351.86757 

-31.0 

117.6 

MP2/6-3 1  G(d)//RHF/6-3  lG(d) 

-352.92725 

-12.2 

106.2 

MP2/6-3  l+G(d)//RHF/6-3  lG(d) 

-352.94539 

-9.8 

108.3 

HNBN3  Imidazidoborane  CSeeEa.  101 

RHF/6-31G(d) 

-243.01229 

-57.6 

40.4 

MP2/6-3 1  G(d)//RHF/6-3  lG(d) 

-243.73651 

-38.0 

60.0 

MP2/6-3  l+G(d)//RHF/6-3  lG(d) 

-243.75290 

-41.1 

56.9 

In  a  fashion  similar  to  B2H6,  TADB  can  dissociate  upon  heating,  yielding  diazidoborane 
(DAB); 

H 

N3  /  \  N3  N3 

\  /  \  /  / 

B— B  ^  2H-B  (9) 

/\  /  \  \ 

N3  \  /  N3  N3 

H 

Ah  initio  calculations  of  diazidoborane  indicate  a  vibrationally  stable  C2v  structure  as  shown  in 
Tables  14  and  15.  The  calculated  energies  in  Table  16  indicate  that  reaction  (9)  is  exothermic  by 
about  9.8  kcal/mol. 

DAB,  like  the  halogen  azides,  can  eliminate  a  terminal  N2  group.  Calculations  at  the 
RHF/6-3  lG(d)  level  of  theory  indicate  that  the  barrier  for  N2  elimination  is  17.9  kcal/mol,  a  value 
slightly  greater  than  that  found  for  FN3  (10.8  kcal/mol  calc.,  13.6  kcal/mol  exp.).  The 

N 

/ 

molecular  fragment  H  -  B  that  results  from  N2  elimination  from  DAB  undergoes  a 

\ 

N3 
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conformational  change  by  migration  of  the  lone  N  atom  away  from  the  N3  group  and  insertion 
into  the  H-N  bond  to  form  a  planar  imidazidoborane  (HNBN3)  structure.  Calculations  at  the 
RHF/6-31G(d)  level  of  theory,  shown  in  Tables  14  and  15,  indicate  that  this  compound  is 
vibrationally  stable.  The  overall  reaction 


N3 

/ 

H  -  B  HNBN3  +  N2 

\ 

N3 

is  exothermic  by  41,1  kcal/mol,  as  calculated  from  the  data  in  Table  16. 


(10) 


3.3  Nitrated  Cyclic  and  Polycyclic  Compounds 

Calculations  were  initiated  on  several  new  cyclic  and  polycyclic  compounds  that  appear  to 
have  potential  in  solid  rocket  propellant  formulations.  These  compounds  are  new  developments 
in  the  chemistry  of  nitramines.  Cyclic  nitramines  such  as  RDX  and  HMX  have  been  known  for 
years  and  were  utilized  as  explosives  during  World  War  H.  Mixtures  of  HMX  and  ammonium 
perchlorate  containing  up  to  70%  HMX  have  been  fired  in  demonstration  motors  at  United 
Technologies  Chemical  Systems  Division  (CSD).  Renewed  interest  in  nitramines  centers  around 
the  development  of  caged  polycyclic  structures.  Starting  from  ring  structures  such  as  RDX, 
HMX  or  a  general  alcane  (Cn,+nH2in+nX  nitramine  caged  structures  can  be  synthesized  by 
forming  C-C  bridges  between  the  rings  of  two  or  more  alcanes.  Thus,  adamantane  (CjoHig)  and 
wurtzitane  (Ci2Hig)  can  be  formed  from  bridging  five-member  and  six-member  rings, 
respectively.  Subsequent  nitration  leads  to  very  energetic  forms  of  these  caged  molecules. 

3.3.1  RDX  and  HMX.  A  study  of  polycyclic  caged  nitramines  was  initiated,  using  RDX 
and  HMX  as  prototype  systems  for  evaluation  and  calibration  of  the  computational  approaches. 
Many  of  the  molecules  investigated  were  too  large  for  applications  of  ai  initio  techniques. 
Therefore,  semi-empirical  methods  such  as  AMI  were  compared  with  ab  initio  Hartree-Fock 
calculations  for  the  smaller  systems  in  order  to  gain  some  confidence  in  the  reliability  of  the  semi- 
empirical  predictions.  The  first  system  examined  was  RDX  in  both  the  boat  and  chair  structures. 
The  optimized  ab  initio  and  AMI  structures  for  the  chair  form  are  compared  in  Figure  8. 

UTRC's  calculated  bond  lengths  and  angles  are  in  good  agreement  with  the  data  compiled  by 
Rogers  (Ref  25).  The  AMI  predictions  are  very  close  to  the  optimized  HF/STO-3G  ab  initio 
structure.  AMI  bond  lengths  are  within  0.05  A  and  AMI  bond  angles  are  within  3%  of  the  ab 
initio  results. 
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The  second  system  examined  was  HMX  where  only  the  crown  (C4y)  structure  was 
examined.  Again  the  AMI  predicted  structure  is  very  close  to  that  obtained  from  a  very  lengthy 
HF/STO-3G  calculation. 


HF/6-31G(d,  p) 
AMI 


Figure  8 

Comparison  of  HF/6-31G(d,  p)  and  AMI  Structure  for  RDX  (chair) 


3.3.2  CL20.  Based  on  observations  that  AMI  predicted  reliable  geometries  for  RDX  and 
HMX,  a  study  was  undertaken  of  polycyclic  caged  structures,  beginning  with  the  six-carbon 
isowurtzitane  structure  [CL20].  In  order  to  further  verify  the  reliability  of  AMI  predicted 
structures,  both  ab  initio  and  AMI  optimized  structures  were  calculated  for  the  zeta  [C2v]  fonn 
of  hexanitrohexaazaisowurtzitane,  a  polycyclic  nitramine  (PCN)  sometimes  referred  to  as  CL20. 
The  ab  initio  calculations  were  carried  out  at  the  HF/STO-3G  level  of  theory,  and  the  AMI 
method  employs  the  STO-3G  basis  for  describing  the  valence  shell.  The  structures  are  compared 
in  Figure  9,  where  some  small  (<5°)  deviations  can  be  seen  in  the  predicted  bond  angles  for  the 
nitro  groups,  but  there  is  overall  excellent  agreement  between  the  ab  initio  and  semi-empirical 
predictions.  This  is  a  36-atom  molecule  and  represents  the  practical  limit  in  the  size  of  the  PCN 
structures  that  were  planned  to  be  investigated.  The  less  symmetric  conformations  of  this  PCN 
were  located  using  the  AMI  method.  They  are  shown  in  Figure  10  for  the  zeta  [C2v]  and  epsilon 
[Cj]  forms.  These  conformations  do  not  necessarily  correspond  one-to-one  with  the  observed 
polymorphs.  In  particular,  it  is  known  that  the  alpha  phase  of  CL20  has  a  non-characteristic 
infrared  (IR)  spectra  indicating  a  solid  rearrangement  or  possibly  a  hydrated  form.  A  companion 
of  UTRC's  calculated  IR  spectra  for  the  C2v  form  of  CL20  with  the  experimental  spectra  reported 
by  Yuxiang,  et  al.  (Ref  26)  is  shown  in  Figures  1 1  and  12,  respectively.  It  is  clear  that  the 
primary  features  corresponding  to  nitro  groups  are  well  represented  by  these  calculations. 
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Figure  9 

Comparison  of  the  Predicted  Structure  of  the  Zeta  Form  of  Isowurtzitane 


zeta  (C2v) 

Figure  10 

Comparison  of  the  Zeta  and  Epsilon  Forms  of  Isowurtzitane  Using  AMI  Theory 
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%  Transmission  Absorption 


Frequency,  (citt'') 


Figure  11 

Calculated  IR  Spectrum  for  PCN/CL20  -  Cjy  Form  (HF/STO-3G  Level  of  Theory) 
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Figure  12 

FTIR  Spectra  of  Four  Polymorphs  of  CL20 


3.3.3  New  PCNs.  Based  on  the  success  of  the  AMI  method  in  describing  CL20,  the 
structure  of  several  new  polycyclic  nitramine  (PCN)  compounds  were  anal5^ed  at  the  AMI  level 
of  theory.  The  molecules  under  study  are  illustrated  in  Figure  13.  The  first  PCN  structure  is 
representative  of  carbon-bridging  on  two  RDX  molecules.  There  is  a  vibrationally  stable 
structure  in  Cg  symmetry.  This  compound  would  have  an  oxygen  balance  of -3  O  atoms.  The 
second  PCN  structure,  shown  in  Figure  13,  is  representative  of  carbon-bridging  on  two  HMX 
molecules.  Again,  there  is  a  vibrationally  stable  structure  in  Cg  symmetry.  The  compound  would 
have  an  oxygen  balance  of -4  O  atoms.  Next,  a  nitramine  group  addition  was  tried  across  the  C-C 
brigeheads  in  CL20.  This  results  in  the  third  compound  illustrated  in  Figure  13,  which  exhibits 
exceptional  vibrational  stability.  This  is  the  first  polycyclic  nitramine  that  has  a  positive  oxygen 
balance  of +3  O  atoms.  Its  thermal  stability  is  unknown,  but  several  synthesis  routes  are  under 
discussion  at  CSD.  A  laboratory  synthesis  program  seems  to  be  warranted  for  this  compound. 

3.4  Dinitramide  Compounds 

Recently,  Bottaro  and  Schmitt^  have  synthesized  several  new  compounds  that  exhibit  a 
dinitramide  anion  structure.  These  compounds  form  as  salts  similar  to  ammonium  nitrate  or 
ammonium  perchlorate,  rather  than  as  the  less  stable  organic  alkyldinitramines  (Ref  27).  UTRC 
has  examined  the  structure  of  this  new  inorganic  anion  and  reported  estimates  of  its  vibrational 
spectra  and  thermodynamic  stability.  The  dinitramide  anion  exhibits  exceptional  stability,  and 
several  new  members  of  the  inorganic  dinitramide  femily  have  been  examined  for  structural 
stability  and  energy  content.  Ab  initio  calculations  were  carried  out  for  several  new  dinitramide 
oxidizers  at  the  HF/6-3  lG(d)  and  MP2/6-3  lG(d)  levels  of  theory. 

3.4.1  DN.  The  structure  of  the  dinitramide  anion  has  been  studied  at  both  the  HF  and  MP2 
levels  of  theory.  The  structure  of  the  anion  is  illustrated  in  Figure  14,  where  calculated  bond 
lengths  and  angles  are  compared  with  data  derived  by  Gilardi^  fi'om  X-ray  analysis  of  ADN.  A 
comparison  of  the  calculated  IR  features  for  the  dinitramide  anion  with  data  taken  on  the  cesium 
salt  is  shown  in  Figure  15.  The  most  prominent  feature  is  the  new  N-N-NO2  rock  motion  at  ~766 
cm"^  which  characterizes  this  species.  The  dinitramide  anion  has  C2  symmetry,  and  this  symmetry 
is  preserved  when  forming  an  alkali  salt  such  as  lithium  dinitramide  (LiDN)  or  cesium  dinitramide 
(CsDN). 


*  Scbmitt,  R.  J.  and  Bottaro,  J.  C.:  private  coimnunication,  SRI  International,  1992. 
^  Gilardi;  private  communication.  Naval  Research  Laboratory,  1992. 
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Figure  13 

New  PCN  Structures 
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(nitrogen  -  bridging  on  2  RDX  molecules) 
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Figure  14 

Calculated  Bond  Lengths  and  Angles  of  the  Dinitramide  Anion  [N(N02)2~] 
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Figure  15 

Characteristic  IR  Features  of  the  Dinitramide  Anion  [N(N02)2*] 
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3.4.2  ADN.  The  structure  of  ADN  was  analyzed  at  both  the  HF  and  MP2  correlated  levels 
of  theory.  In  Figure  16  the  optimum  structure  for  ADN  at  the  RHF/6-3  lG(d,  p)  level  is  shown. 
This  structure  is  very  similar  to  that  found  by  Gilardi^  for  the  solid  material  and  comprises  a 
NH4'^(N02)2’  salt-like  structure.  A  second  stable  form  of  ADN,  as  illustrated  in  Figure  17,  is 
found  at  both  the  RHF  and  MP2  levels  of  theory.  This  structure,  which  probably  represents  the 
gas  phase  form  of  this  material,  consists  of  a  weak  NH3-N(N02)(N00H)  acid-base  complex. 

The  decomposition  profiles  observed  by  Fetheroff  and  Litzinger^  and  by  Snelson^  are  consistent 
with  decomposition  of  this  acid-base  complex  to  yield  NH3  and,  through  a  concerted 
unimolecular  decomposition  of  the  N(N02)2  group,  to  yield  N2O,  H2O,  N2  and  NO2. 
Thermodynamic  calculations  of  ADN  are  incomplete,  but  preliminary  calculations  indicate 
AHf(s)(298)  =  -36  kcal/mol  for  ADN,  in  agreement  with  experimental  data  from  China  Lake.^ 


ADN  [NH4N(N02)2l  -  Cj  Structure  -  RHF/6-31G(d,  p)  Level  of  Theory 


ADN  INH3-HN(N02)2  Complex]  -  Cg  Structure  -  MP2/6-31G(d)  Level  of  Theory 


^  Gilardi:  private  communication.  Naval  Research  Laboratory,  1992. 

^  Fetheroff,  B.  L.  and  Litzinger,  T.  A.;  private  communication,  Pennsylvania  State  University,  1992. 

4 

Snelson,  A.:  private  communication,  IIT  Research  InsUtute,  1993. 

*  Nadler;  private  communication.  Naval  Weapons  Center,  China  Lake,  CA,  1995. 
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3.4.3  HADN.  The  improved  oxygen  balance  ofhydroxyammonium  nitrate  (HAN)  over 
ammonium  nitrate  suggested  an  examination  of  the  stability  of  the  parallel  HADN  structure.  It 
was  found  that  HADN  is  stable  as  a  C|  ionic  structure  with  a  calculated  heat  of  formation  of 
~  -30  kcal/mol.  Recent  s)mthesis  and  characterization  studies  of  this  compound  at  SRI  indicate 
that  it  is  a  liquid  under  normal  conditions  of  temperature  and  pressure.  Further  studies  of  HADN 
were  also  carried  out  at  the  MP2  level  of  theory  to  more  fiilly  optimize  the  molecular  structure. 
The  results  are  illustrated  in  Figure  18,  which  suggests  a  salt-like  structure,  similar  to  that 
observed  for  ADN.  Preliminary  laboratory  data  from  SRI  suggest  that  HADN  hydrolyzes  readily 
and  is  probably  only  stable  at  low  temperatures. 


Figure  18 

HADN  [NH3OH  N(N02)2l  -  Ci  Structure  -  RHF/6-31G(d,  p)  Level  of  Theory 


3.4.4  Dinitraminic  Acid.  In  the  dinitramide  series,  of  particular  interest  is  the  simplest 
compound,  hydrogen  dinitramide,  which  could  exist  as  the  smallest  member  of  an  organic  alkyl- 
dinitramine  series,  R-N(N02)2,  or  possibly  as  an  inorganic  acid  of  the  form,  N(N02)(NOOH). 
Leroy,  et  al.  (Refs.  28,  29)  have  reported  gradient  optimized  RHF/4-3 IG  and  RHF/6-3 IG 
energies  for  HN(N02)2  and  have  made  semi-empirical  estimates  of  its  heat  of  formation,  but  they 
have  not  characterized  the  calculated  molecular  conformation  or  vibrational  stability. 

Ab  initio  calculations  on  HN(N02)2  were  performed  at  both  the  HF  and  MP2  levels  of 
theory  using  the  standard  6-3  lG(d,  p),  6-3 1+G(d,  p)  and  6-3 1 1+G(d,  p)  basis  sets.  Although  the 
calculated  geometries  are  quite  insensitive  to  the  addition  of  diffuse  basis  functions,  their  use  is 
necessaiy  for  thermochemical  calculations  involving  the  anion.  Using  the  Gaussian  system  of 
programs  (Ref  15),  gradient  optimizations  were  performed  to  locate  the  geometries  of  the 
stationary  points.  Harmonic  vibrational  frequencies  were  then  computed  to  verify  that  the 
calculated  structures  are  true  minima.  The  CADPAC  electronic  structure  code  (Ref  30)  was 
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used  to  perform  the  frequency  analysis,  since  analytic  second  derivatives  are  available  at  the  MP2 
level  of  theory.  Two  stable  conformations  were  located  at  both  levels  of  theory.  To  assist  in  the 
analysis  of  the  gas  phase  acidity  of  this  molecule,  parallel  level  of  theory  calculations  were  carried 
out  for  the  dinitramide  anion.  This  anion  is  found  to  be  a  C2  structure  at  all  levels  of  theory.  The 
results  are  summarized  in  Appendix  B-3,  which  contains  an  abstract  of  a  technical  paper  on  this 
system. 

The  minimum  energy  stationary  point  on  the  ground  state  surface  corresponds  to  a  Cg 
structure  of  the  secondary  amine  form  (nitro-form)  [HN(N02)2]  of  hydrogen  dinitramide.  A 
calculation  of  the  4-3 IG  optimized  geometry  indicates  that  this  is  the  structure  found  by  Leroy,  et 
al.  (Refs.  28, 29).  A  second  planar  Cg  structure  is  a  metastable  acid  isomer.  This  acid-form 
[N(N02)(N00H)]  of  hydrogen  dinitramide  is  found  to  lie  9.4  kcal  higher  in  energy  than  the  nitro- 
form.  Two  other  stable  acid-forms  are  found  which  correspond  to  non-planar  Cg  structures 
arising  from  rotation  of  the  nitro  group.  The  chemistry  of  hydrogen  dinitramide  thus  appears  to 
be  similar  to  that  observed  for  nitromethane,  [CH3NO2],  where  two  tautomeric  forms  are  known 
to  exist  (Ref  3 1).  Deprotonation  of  nitromethane  occurs  slowly  in  basic  solution,  and  upon 
reacidification,  the  strongly  acidic  isomer  is  formed.  In  both  compounds  the  acid-form  reverts  to 
the  more  stable  nitro-form.  The  time  constant  for  conversion  between  the  two  forms  is  dependent 
on  the  activation  energy  of  the  transition  state  above  the  acid-form.  Apparently  there  is  a  sizable 
barrier  for  nitromethane,  since  the  conversion  occurs  slowly  at  room  temperature. 

3.4.5  Trinitramine.  Although  there  are  many  possible  compounds  of  the  form  RR'NN02 
and  RN(N02)2,  there  is  only  one  possible  trinitramine  [N(N02)3].  There  has  been  speculation 
concerning  the  existence  of  trinitramine,  but  very  few  results  have  actually  been  reported  to  date. 

Ab  initio  calculations  were  performed  on  pyramidal  (C3)  and  planar  (D3ji)  structures  of 
N(N02)3.  Full  geometry  optimizations  and  harmonic  vibrational  frequency  calculations  were 
performed  on  both  planar  and  pyramidal  structures  at  the  RHF/6-3  lG(d)  level  of  theory.  An 
additional  MP2/6-31G(d)  geometry  optimization  and  harmonic  vibrational  frequency  calculation 
was  performed  on  the  pyramidal  structure.  Thermochemical  calculations  include  vibrational  zero- 
point  effects,  estimated  from  the  RHF/6-3  lG(d)  vibrational  frequencies  scaled  by  0.8929.  A 
summary  of  the  calculated  results  is  given  in  Tables  17  and  1 8.  The  C3  structure  is  shown  in 
Figure  19. 
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Table  17.  Calculated  Energies  and  Geometries  of  Trinitramine  [N(N02)3] 

Energies  are  in  hartrees  and  geometries  are  in  angstroms  and  degrees 


Level  of  Theory 

RHF/6-31Gfd) 

MP2/6-31Gld) 

RHF/6-31Gfd') 

Symmetry 

C3 

C3 

^3h 

Enersv^ 

-66.44582  (0) 

-668.26537  (0) 

-66.40519 

Geometry 

N-N 

1.4444 

1.5453 

1.4663 

N-Oi 

1.1713 

1.2191 

1.1698 

N-O2 

1.1757 

1.2209 

1.1698 

N-N-N 

110.9 

105.2 

120.0 

N-N-Oj 

112.6 

111.0 

116.2 

N-N-O2 

116.7 

117.4 

116.2 

N-N-N-Oj 

-70.8 

-76.0 

0.0 

N-N-N-O2 

110.4 

105.2 

180.0 

®  The  number  of  imaginary  vibrational  frequencies  is  indicated  in  parenthesis. 


Table  18.  Calculated  Vibrational  Frequencies  and 
IR  Intensities  of  Trinitramine  [N(N02)3] 


Level  of  Theory  RHF/6-31G(d) 

Symmetry  C3 


MP2/6-31G('d) 

C3 


Frequencies  are  in  cm'^  and  IR  intensities  [  ]  are  in  km/mol 


Group 


A  E 


127 

[0.6] 

53 

[0.1] 

241 

[5.0] 

249 

[3.7] 

444 

[0.0] 

508 

[2.3] 

651 

[4.7] 

801 

[8.7] 

937 

[10.3] 

908 

[155.5] 

997 

[9.3] 

1155 

[79.3] 

1641 

[31.5] 

1537 

[290.8] 

1934 

[869.1] 

1979 

[550.1] 

A 

E 

96 

[0.1] 

43 

[0.1] 

217 

[2.1] 

197 

[4.2] 

325 

[0.1] 

321 

[58.0] 

521 

[0.2] 

605 

[59.7] 

796 

[12.1] 

748 

[172.1] 

835 

[6.0] 

800 

[70.9] 

1350 

[14.3] 

1263 

[180.6] 

1912 

[90.5] 

1921 

[58.9] 

32 


o 


Figure  19 

Trinitramine  Molecule  [N(N02)3]  -  C3  Structure  -  MP2/6-31G(d)  Level  of  Theory 


As  in  anunonia,  nitrogen  trifluoride,  and  trimethylene,  calculations  show  that  trinitramine  is 
pyramidal,  not  planar.  At  the  RHF/6-31G(d)  level,  the  planar  structure  is  24  kcal/mol  higher  in 
energy  than  the  pyramidal  ground  state.  The  calculated  vibrational  frequencies  for  the  pyramidal 
structure  are  all  real,  indicating  this  structure  is  a  true  minimum.  Three  of  the  RHF/6-31G(d) 
vibrational  frequencies  found  for  the  planar  trinitramine  stracture  are  imaginary,  indicating  that  it 
is  not  a  stable  minimum  or  a  transition  state,  and  therefore  no  further  calculations  were 
performed  on  this  structure. 

The  heat  of  formation  of  trinitramine,  can  be  estimated  from  the  heat  of  reaction  for 

NH3  +  3HNO3  ^  N(N02)3  +  3  H2O,  AHf®  = +63.1  kcal/mol  (11) 

or,  alternatively,  from 

3NH2NO2  ^  N(N02)3  +  2NH3,  AHf®  = +55.0  kcal/mol  (12) 

The  MP2/6-31G(d)  heats  of  formation  calculated  from  the  two  reactions  were  averaged  to  obtain 
the  final  estimate,  AHf°(N405)  =  +59.0  kcal/mol.  Analogous  calculations  (Ref.  20)  of  the  heat  of 
formation  of  nitric  acid  give  an  error  of  1.7  kcal/mol.  UTRC's  estimated  error  in  the  calculated 
heat  of  formation  of  trinitramine  is  ±5  kcal/mol. 
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3.4.6  Nitrohydrazine  Compounds.  The  nitration  of  hydrazine  has  been  examined  since 
the  chemistry  parallels  that  of  the  nitramides.  Three  vibrationally  stable  compounds  are  found 
from  systematic  nitration  of  hydrazine:  1,1-dinitrohydrazine  [aminodinitramide],  1,2- 
dinitrohydrazine  and  tetranitrohydrazine  [dinitramide  dimer].  The  first  and  last  of  these 
compounds  are  clearly  in  the  dinitramide  family.  1,2-dinitrohydrazine  exhibits  the  highest  stability 
in  this  group,  with  a  central  N-N  bond  length  of  1 .36  A.  All  of  these  compounds  are  highly 
energetic  and  are  suggested  target  molecules  for  synthesis  activities.  The  structure  of  these 
molecules  is  illustrated  in  Figure  20  (see  page  35).  Tetranitrohydrazine  [(N(N02)2)2]  with  an 
estimated  heat  of  formation  of +128  kcal/mol  and  typical  composite  propellant  density  equal  to 
that  for  AP,  remains  the  main  target  compound  for  synthesis  in  this  group. 

3.4.7  Dinitroborane.  Studies  were  carried  out  on  several  nitroborane  structures,  including 
an  extensive  analysis  of  the  simplest  such  compound,  BH2NO2.  Six  structures  of  BH2NO2  have 
been  theoretically  characterized.  The  calculated  results  are  summarized  in  Table  19  and 

Figure  21 .  Structures  1  and  2  have  been  previously  studied  (Ref  32),  while  Structures  3-6  are 
new. 


Table  19.  Calculated  Results  for  BH2NO2  Structures 


SUucture 

RHF/6-3  IGfd) 

MP2/6-31Gfdi 

1 

-229.87736  [IJ 

-230.49038  [1] 

2 

-229.88688  (0] 

-230.49858  [1] 

3 

— 

-230.50913  [0] 

4 

-229.92361  [1] 

-230.53460  [0] 

5 

-229.93967  [0] 

-230.53965  [0] 

6 

-229.93285  [1] 

-230.53598  [0] 

The  number  of  imaginary  vibrational  frequencies  is  given  in  parenthesis  following  the  calculated  energy.  No  RHF 
stationary  point  was  found  for  Structure  3. 


From  Figure  21,  Structures  1-3  can  be  seen  to  have  a  nitro-form,  and  Structures  4-6  are  nitrite 
compounds.  Electron  correlation  plays  a  very  significant  role  in  determining  the  vibrational 
stability  of  the  various  BH2NO2  structures,  as  three  of  the  six  reported  structures  change  from 
minima  to  saddles,  or  vice  versa,  when  the  level  of  theory  is  taken  from  RHF/6-3  lG(d)  to  MP2/6- 
3  lG(d).  This  is  a  peculiarity  of  nitroboranes,  because  of  their  electron  deficient  nature.  Also  a 
consequence  of  the  electron  deficiency  is  the  tendency  to  form  bridged  Structures  (3-4)  and  to 
prefer  the  formation  of  B-0  bonds.  Thus,  the  stable  structures  are  somewhat  different  than  those 
found  for  nitromethane  and  nitramide. 
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1,2-dinitrohydrazine 

['A]  -  C2  Symmetry 


Tetranitrohydrazine 

[*A]  -  C2  Symmetry 


Figure  20 

Nitrohydrazine  Compounds 
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Figure  21 

BH2NO2  Structures 
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Heats  of  formation  can  be  estimated  from  the  reaction 


BH3+HNO3  H2O  +  BH2NO2  (13) 

using  the  calculated  MP2/6-31G(d)  heat  of  reaction  and  the  known  heats  of  formation  ofborane, 
nitric  acid,  and  water.  The  heat  of  formation  of  Structure  3  is  estimated  to  be  -4  ±5  kcal/mol  and 
that  of  Structure  5  is  estimated  to  be  -16  ±5  kcal/mol. 

A  transition  state  for  the  interconversion  of  Structures  3  and  4  has  been  located  and  gives  a 
barrier  >100  kcal/mol  for  this  process.  Therefore,  it  may  be  that  both  nitro  and  nitrite  forms  of 
BH2NO2  exist  as  stable  species.  The  interconversion  of  Structures  4-6  is  undoubtedly  more 
facile. 

Further  studies  at  higher  levels  of  theory  will  be  necessary  to  determine  conclusively  the 
energetics  and  stability  of  BH2NO2.  However,  it  appears  from  the  present  calculations  that  there 
are  no  direct  boron  analogues  of  nitromethane  and  nitramide,  and  this  suggests  that  boron 
analogues  of  dinitramide  will  probably  be  bridged  or  other  such  unconventional  structures.  This 
speculation  is  confirmed  by  detailed  calculations  on  dinitroborane  structures.  The  HB(N02)2 
analogue  of  HDN  is  found  to  be  a  saddle  point  at  the  MP2/6-3  lG(d)  level  of  theory,  but  a  bridged 
C|  structure,  analogous  to  BH2NO2,  Structure  3,  is  found  to  be  a  stable  minimum.  This  structure 
is  shown  in  Figure  22.  Additionally,  a  C2v  ring  structure,  shown  in  Figure  23,  has  been  found  to 
be  stable  at  the  MP2/6-3  lG(d)  level,  and  lies  about  60  kcal/mol  higher  in  energy  than  the 
structure  of  Figure  22. 

Further  work  is  needed  to  characterize  these  structures,  and  assess  their  potential  as  high 
energy  density  matter  (HEDM)  compounds. 
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Figure  23 

[HB(N02)2]  -  C2v  Structure 


3.5  Nitrated  Ring  Structures 

3.5.1  NTO.  A  study  was  initiated  on  substituted  triazole  and  tetrazole  compounds  that 
may  exhibit  high  energy  content  with  less  sensitivity  to  shock  than  nitrated  cyclic  compoimds 
such  as  RDX  and  HMX.  The  development  of  explosives  with  lessened  sensitivity  to  detonation 
by  shock  or  impact  is  an  area  of  considerable  current  interest.  There  are  several  factors  which, 
when  taken  together,  determine  the  impact  or  shock  sensitivity  of  an  energetic  material.  At  one 
end  of  the  scale  is  the  interaction  of  the  bulk  material  (crystal  or  polymer)  with  the  applied 
external  force,  leading  to  deformation,  compression,  and  ultimately,  the  localization  of  energy, 
creating  "hot  spots".  At  the  other  end  of  the  scale  is  the  molecular  thermochemistry  of  the 
energetic  material,  as  measured  by  bond  dissociation  energies  and  barrier  heights  for 
rearrangement  and  decomposition.  These  determine  the  molecular  pathways  for  assimilation  of 
energy  applied  externally  to  the  bulk  material.  Full  understanding  of  the  natme  of  shock 
sensitivity  in  energetic  materials  requires  investigation  of  both  bulk  and  molecular  phenomena. 

In  particular,  3-nitro-l,2,4-triazol-5-one  (NTO)  is  now  being  considered  for  use  as  an 
insensitive  explosive.  Partly  based  on  semiempirical  (PM3)  calculations,  preliminary  estimates 
have  been  made  of  the  thermochemistry  of  the  thermal  decomposition  of  NTO.  The  possible 
pathways  that  have  been  examined  are: 


1)  C-NO2  bond  breaking 


AH  ~  65+10  kcal/mol 
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2)  NO2  ONO  rearrangement  and  subsequent  bond  breaking 
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AH^  ~  75+10  kcal/mol 


AH  ~  78+10  kcal/mol 


4)  concerted  unimolecular  decomposition 


O 

X 

HN^^NH 

N=^ 

NO2 


NH  +  HNCO  +  NCNO, 


5)  nitronic  acid  formation 


AH^  ~100±10  kcal/mol 
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A. 


HN  NH 

K 

NO2 


OH 


AH^  >100  kcal/mol 


UTRC's  estimate  of  65  kcal/mol  for  the  C-NO2  bond  is  about  20-30  kcal/mol  above  current 
estimates  of  the  N-NO2  bond  strength  in  materials  such  as  RDX  and  provides,  in  combination 


39 


with  thermochemical  estimates  for  Alternatives  2  through  5,  a  simple  molecular  explanation  of  the 
lessened  sensitivity  of  NTO  relative  to  RDX  and  HMX. 

Although  the  lowest  energy  pathway  suggested  by  preliminary  calculations  is  via  Alternative 
1  (C-NO2  bond  breaking).  Alternatives  2  and  3  are  energetically  close  enough  that  more  accurate 
calculations  might  reverse  their  order.  Indeed,  there  is  evidence  that  NTO/acetone  solutions 
decompose  via  solvent-assisted  hydrogen  abstraction,  and  it  is  possible  that  N-H  bond  breaking 
could  be  the  lowest  energy  decomposition  channel.  It  is  clear  that  a  more  detailed  study  of  the 
molecular  thermochemistry  of  NTO  and  similar  materials  is  necessary  to  understand  their  thermal 
decomposition  pathways. 

3.5.2  Nitrozoles.  Studies  were  initiated  on  the  structure  and  stability  of  several  substituted 
triazoles  and  tetrazoles  based  on  discussions  with  the  SRI  synthesis  group.  The  baseline 
compound  for  comparison  is  NTO,  which  is  a  shock-insensitive  explosive.  NTO  is  described  in 
detail  in  the  previous  section.  There  is  a  stable  structure  with  Cg  symmetry  for  1 -amino-3 -nitro- 
1,2, 4-triazole  (ANT);  1 -amino-3, 5 -dinitro- 1,2,4-triazole  (ADNT);  and  a  planar  structure  with  Cg 
symmetry  for  l-fluoro-3-nitro- 1,2,4-triazole.  All  of  these  compounds  are  highly  energetic  relative 
to  NTO. 

Higher  energy  structures  can  be  synthesized  by  further  nitrogen  substitution  in  the  ring.  For 
example,  there  is  stability  for  1 -amino-3 -nitrotetrazole  in  Cg  symmetry.  This  compound  has  a  very 
high  heat  of  formation  (+96.5  kcal/mol)  and  has  recently  been  synthesized  at  SRI. 

The  calculated  structures  for  this  class  of  energetic  compounds  are  illustrated  in 
Figures  24  through  28.  The  calculated  solid  densities  and  estimated  gas  phase  heats  of  formation 
are  shown  in  Table  20. 


Figure  24 

[NTO]  3-nitro-l,2,4-triazole-5-one  (C2H2N4O3)  -  Cg  Symmetry,  *A’ 
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Figure  25 

[ANT]  l-amino-3-iiitro-l,2,4-triazoIe  (C2H3N5O2)  - 


[ADNT]  l-aniino-3,5-dinitro-l^, 4-triazole  (C2H2N6O4)  - 


Cg  Symmetry,  *A* 


Cj  Symmetry,  ^A’ 
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Figure  27 

[FNT]  l-fluoro-3-nitro-l, 2, 4-triazole  (C2HFN4O2)  -  C*  Symmetry,  ^A' 


Figure  28 

l-amino-3-nitrotetrazole  (CH2N6O2)  -  Cg  Symmetry,  *A* 


Table  20.  Thermochemistry  of  Advanced  Propellants 


Density 

AHf(25°C) 

Comoound 

Formula 

Mol.  Wt. 

em/cm^ 

kcal/mol 

1.80 

-36.0 

ammonium  dinitramide  [ADN] 

H4N4O4 

124.056 

1.75 

+8.3  (g) 

nitrosyldinitramide 

N4O5 

136.023 

1.94 

+49.8 

trinitramine 

N4O6 

152.022 

2.15 

+80.4 

1,1-dinitrohydrazine  (aminodinitramide) 

H2N4O4 

122.040 

2.04 

+58.6 

1 ,2-dinitrohydrazine 

H2N4O4 

122.040 

2.04 

+44.3 

tetranitrohydrazine 

N608 

212.034 

2.17 

+128.2 

bis  [nitroamino]  nitrosylamine  anion 

N505- 

150.030 

1.92 

+57.9 

diazidamine 

H1N7 

99.057 

1.56 

+199.0 

triazidamine 

Nio 

140.070 

1.58 

+302.0 

diazidonitronium  cation 

N70'^ 

114.048 

1.65 

+397.1 

nitrosylazide 

N40 

72.027 

1.54 

+115.4 

nitrylazide 

N402 

88.026 

1.76 

+95.1 

3-nitro-l,2,4-triazol-5-one  [NTO] 

C2H2N403 

130.063 

1.75 

0 

1 

l-amino-3-nitro- 1,2, 4-triazole  [ANT] 

C2H3N502 

129.079 

1.61 

+55.6 

[1.70] 

[+50.0]  (exp.) 

1  -amino-3 ,5-dinitro- 1 ,2,4-triazole  [ADNT] 

C2H2N604 

174.076 

1.77 

+73.2 

5-amino-2,4-dinitro- 1 ,3,5-triaza- 

C2H2N604 

174.076 

1.92 

+173.3 

bicyclo[2, 1 ,0]pent-2-ene 

l-fluoro-3-nitro- 1,2,4-triazole  [FNT] 

C2HFN402 

132.054 

1.68 

+76.2 

1  -amino-3-nitrotetrazole 

CH2N602 

130.067 

1.67 

+96.5 

dioxopentazolide  anion 

N502- 

102.033 

1.76 

+100.4 

nitroacetylene 

C2HN02 

71.035 

1.51 

+71.7 

dinitroacetylene 

C2N204 

116.032 

1.79 

+94.8 

1-nitropropylene 

C3H3N02 

85.062 

1.37 

+57.8 

nitrylcyanide 

CN202 

72.023 

1.82 

+58.0 

1 ,3,5-trinitro- 1 ,3 ,5-triazacyclohexane  [RDX] 

C3H6N606 

222.117 

1.76 

+53.5 

1 ,3-dinitro- 1 ,3-diazacyclobutane  [HEDMX] 

C2H4N404 

148.078 

1.76 

+53.7 

1,3,3  -trinitroazididine  [TN AZ] 

C3H4N406 

192.087 

1.78 

+35.0 

43 


3.5.3  Nitramines/Azididines.  Cyclic  nitramines  that  also  contain  geminal  dinitroalkyl 
groups  have  recently  been  synthesized  (Refs.  33-36).  These  compounds  may  exhibit  an  increase 
in  oxygen  balance  over  the  corresponding  compounds  with  mononitroalkyl  groups.  A  compoimd 
of  primary  interest  in  this  group  is  TNAZ,  a  four-membered  ring  in  the  RDX,  HMX  series  that 
exhibits  a  single  dinitroalkyl  group.  Other  compounds  in  this  series  include  1, 3,5,5- 
tetranitrohexahydropyrimidine  [DNNC]  and  1,3,3, 5,7, 7-hexanitro-l,5-diazacyclooctane 
[HNDC],  which  are  analogues  of  RDX  and  HMX,  respectively. 

An  extensive  set  of  ab  initio  calculations  has  been  carried  out  on  TNAZ.  Eight  possible 
conformations  of  TNAZ  were  examined  at  the  HF/3-21G  and  HF/6-31G(d)  levels  of  theory.  The 
ground  state  conformation  has  the  nitro  groups  at  ring  Positions  1  and  3  in  a  perpendicular  trans¬ 
arrangement,  relative  to  the  azididine  ring.  A  second  very  low-lying  conformation  (+0.3  kcal) 
has  the  nitro  groups  at  ring  Positions  1  and  3  in  a  perpendicular  cis-arrangement  to  the  azididine 
ring.  There  appears  to  be  nearly  free  rotation  of  the  nitro-groups  at  ring  Position  3.  A  transition 
state  energy  is  calculated  at  +2.0  kcal  when  both  nitro  groups  at  Position  3  are  perpendicular  to 
the  azididine  ring.  Figures  29  and  30,  respectively,  show  these  low-lying  cis-  and  trans¬ 
conformations  of  azididine.  If  similar  results  are  found  at  a  correlated  level  of  theory,  there  may 
be  more  than  one  polymorph  of  this  compound. 


(nitro  groups  at  ring  Positions  1  and  3  are  cis  to  ring) 

Figure  29 

1,3,3-triiiitroazididine  (cis)  -  Cj  Symmetry,  ^A' 
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(nitro  groups  at  ring  Positions  1  and  3  are  trans  to  ring) 

Figure  30 

l^yS-trinitroazididine  (trans)  -  Cg  Symmetry,  ^A* 


Correlated  calculations  of  the  ground  state  of  TNAZ  were  next  attempted.  This  system, 
however,  proved  to  be  too  extensive  for  current  computer  capabilities.  Instead,  based  on  a  fully 
optimized  Cg  structure  (25  geometric  parameters),  the  heat  of  formation  of  this  molecule  was 
calculated  using  the  method  outlined  by  Ibrahim  and  Schleyer  (Ref  37).  A  value  of  AHf^  (25  °C) 
=  +  35.0  kcal/mol  was  obtained  for  the  gas  phase  heat  of  formation.  Estimated  corrections  for 
the  sublimation  energy  indicate  that  TNAZ  is  more  energetic  than  NTO  but  less  energetic  than 
the  recently  synthesized  l-amino-3-nitro- 1,2,4  triazole  compound.^  UTRC's  calculated  density 
for  TNAZ  is  1 .78  g/cc  (Cady's  method),  slightly  greater  than  that  for  NTO  (1 .75  g/cc). 

In  addition  to  studies  of  TNAZ,  an  evaluation  of  cyclodimethylenedinitramine  (1,3-dinitro- 
1 ,3-diazacyclobutane)  was  undertaken.  This  compound  is  the  smallest  ring  member  in  the  RDX, 
HMX  nitramine  series.  There  is  a  stable  structure  with  Cg  symmetry  as  illustrated  in  Figure  3 1 . 
The  estimated  heat  of  formation  of  this  new  compormd  is  +54  kcal/mol,  indicating  that,  on  a 
mass  basis,  this  is  a  much  more  energetic  structure  than  either  RDX  or  HMX.  The  calculated 
solid  density  is  1 .76  g/cc  (Cady's  method),  a  value  somewhat  smaller  than  that  for  RDX  (1 .80 
g/cc). 

Correlated  calculations  of  cyclodimethylenedinitramine  were  undertaken  for  this  smallest 
ring  member  in  the  RDX,  HMX  nitramine  series.  A  stable  C2v  structure  with  *Aj  symmetry  is 
found  at  the  MP2/6-31G(d)  level  of  theory.  The  lowest  frequency  mode  [aj  =  64  cm’^]  is  a 
symmetric  axial  stretch  of  the  C-N-C-N  ring,  indicating  that  this  compound  may  have  lower  ring 
strain  energy  than  expected,  offering  some  chance  for  ring  closure  in  synthesis  attempts.  The 
name  HEDMX  is  suggested  as  an  identifier  for  this  new  cyclic  nitramine. 


*  Schmitt, R.:  private  communication,  SRI  International,  1995. 
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(smallest  member  of  RDX,  HMX  nitramine  series) 
Figure  31 

l,3-dinitro-l,3-diazacycIobutane  -  C2v  Symmetry,  ^ Ai 


3.6  Substituted  Prismanes 

Ab  initio  calculations  were  performed  to  determine  the  structure  and  stability  of 
[3]-prismane  and  its  nitro  derivatives.  (The  simplifying  notation  of  prismane  is  adopted  for  the 
smallest  member  of  the  [n]-  prismane  group  of  compounds.)  In  Figure  32,  the  stable  03^ 
structure  of  unsubstituted  prismane  is  illustrated.  Next  the  prismane  carboxylic  acid  structure 


Figure  32 

Prismane  [C5H5]  as  a  D3],  Structure 
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was  examined,  since  such  substituted  compounds  are  known  for  [4]-prismane  (cubane)  and  may 
provide  a  synthesis  entry  point  for  other  functional  groups.  A  stable  dicarboxylic  acid  structure 
of  prismane  has  been  found  that  exhibits  C2  symmetry,  as  illustrated  in  Figure  33.  In  addition, 
stable  forms  of  both  nitroprismane  and  dinitroprismane  were  found,  as  illustrated  in  Figures  34 
and  35,  respectively. 


Figure  33 

Prismane  Dicarboxylic  Acid  [C5H4(COOH)2l  -  C2  Structure 


Nitroprismane  [C5H5NO2]  -  Q  Structure 
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Figure  35 

Dinitroprismane  [C6H4(N02)2]  -  C2  Structure 

A  preliminary  evaluation  of  the  thermochemistry  of  prismane  and  substituted  prismane  has 
been  carried  out.  Gradient  optimized  structures  for  prismane,  nitrotriprismane,  and 
dinitroprismane  were  obtained  at  the  RHF/6-3  lG(d)  level  of  theory.  Calculations  of  the 
harmonic  vibrational  frequencies  demonstrated  these  to  be  vibrationally  stable  structures. 
Following  Schulman  and  Disch  (Ref.  38),  the  homodesmic  reactions 


6  isobutane 

-> 

prismane  +  9  ethane 

(14) 

5  isobutane  +  nitrobutane 

nitroprismane  +  9  ethane 

(15) 

4  isobutane  +  2  nitrobutane 

dinitroprismane  +9  ethane 

(16) 

were  used  to  calculate  the  heat  of  formation  of  these  structures.  The  enthalpy  change  in  these 
reactions  provides  a  theoretical  estimate  of  the  strain  energy  (Ref  39).  These  calculations 
suggest  that  nitroprismane  and  dinitroprismane  are  not  significantly  more  strained  than  prismane. 
The  calculated  heats  of  formation  and  strain  energies  are  shown  in  Table  21 . 

Table  21.  Calculated  Heats  of  Formation  and  Strain  Energy  of  Prismanes 

MP2/6-31G(d,  p)  Level  of  Theory 


Heat  of  Formation  Strain  Energy 

tkcal/mol')  (kcal/mol) 


Prismane  +135  +147 

Nitroprismane  + 125  +150 

Dinitroprismane  + 122  +159 


Several  synthesis  routes  for  the  nitroprismanes  are  being  evaluated  at  United  Technologies 
Chemical  Systems  Division.  A  five-step  synthesis  route  to  the  dicarboxylic  acid  structure  shown 
in  Figure  33  has  been  formulated.  Formation  of  the  nitro-forms  from  the  dicarboxylic  acid 
structure  may  follow  along  routes  similar  to  those  known  to  jdeld  nitrocubanes  starting  fi'om  the 
carboxylic  acid  form  of  cubane. 


4.0  CONCLUSIONS  AND  RECOMMENDATIONS 

The  main  objective  of  this  research  program  was  to  identify  those  light  element  molecular 
species  that  have  sufficient  internal  energy  to  be  useful  in  advanced  chemical  propulsion  systems. 
Other  objectives  were  to  characterize  the  available  energy  content  of  these  new  compounds,  to 
investigate  probable  kinetic  routes  for  their  formation  and  decomposition,  to  characterize 
vibrational  spectra  and  other  molecular  properties  and  to  suggest  possible  synthesis  routes.  The 
utility  of  these  new  compounds  in  chemical  propulsion  systems  was  anal)^ed  through  calculations 
of  specific  impulse  of  propellant  formulations  that  incorporated  these  materials. 

In  a  previous  study  (Ref.  40)  the  focus  was  on  light  element,  high  symmetry  structures  such  as 
Li3H,  azide-like  structures  such  as  a-N202  and  FN3,  and  hypervalent  structures  such  as  NF5. 
During  the  course  of  the  current  program,  the  emphasis  shifted  to  studies  of  advanced  oxidizers, 
including  characterization  of  the  new  family  of  dinitramide  compounds.  Toward  the  end  of  this 
contract,  as  computational  tools  improved,  it  was  possible  to  study  some  very  complicated 
energetic  ring  structures.  The  compounds  that  were  studied  can  be  categorized  into  six  groups: 

1.  high  symmetry  hypervalent  structures  [OF5] 

2.  azide-like  structures  [N(N3)3,  HN(N3)2,  N(N3)2“,  N(N3)4''',  NON3,  NO2N3,  B2H2(N3)4, 
HB(N3)2,  HNBN3] 

3.  nitrated  cyclic  and  polycyclic  nitramine  compounds  [RDX,  KMX,  CL20,  PCNs] 

4.  dinitramide  compounds  [N(N02)2",  LiDN,  ADN,  HADN,  HN(N02)2,  N(N02)3,  N(N02)4^, 

H2NN(N02)2,  (HNN02)2,  (N(N02)2)2,  BH2NO2,  B(N02)2;  HB(N02)2] 

5.  nitrated  ring  compounds  [NTO,  TNAZ,  ANT,  ADNT,  FNT,  HEDMX] 

6.  substituted  prismanes  [triprismane,  nitrotriprismane,  dinitrotriprismane,  triprismane  carboxylic 
acid] 

A  summary  of  the  results  and  recommendations  for  further  research  into  these  several  groups  of 
compounds  follows. 
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4.1  High  Symmetry  Hypervalent  Structures 

Studies  of  OFg  indicate  that  this  high  symmetry  structure  may  exist  as  a  gas  phase  species 
exhibiting  mainly  covalent  bonding.  However,  there  is  a  low  transition  state  barrier  for  formation 
of  the  ionic  F‘OF5^  structure  which  has  low  energy  content.  The  chemistry  of  this  compound 
thus  resembles  that  for  NF5.  It  may  be  of  interest  to  study  OF4  and  other  mixed  oxidizers,  but  the 
potential  for  significant  improvement  from  the  class  of  oxidizers  appears  to  be  limited. 

4.2  Azide-like  Structures 

All  of  the  new  azide-like  structures  that  were  examined  exhibit  very  large  positive  heats  of 
formation.  The  chemistry  of  the  diazidamide  anion  parallels  that  for  the  dinitramide  anion  but 
exhibits  much  higher  energy  content.  The  fundamental  problem  with  azidammes  is  their  low 
energy  barrier  for  nitrogen  elimination.  This  characteristic  was  observed  experimentally  in  the 
work  by  Benard  at  Rockwell  Science  Center  (Ref  41),  where  TADB  dissociated  upon  heating  to 
yield  B-N  structures  with  much  lower  energy  content.  Nitrosyl  azide  and  nitryl  azide  are  known 
energetic  compounds  but  require  very  low  temperature  synthesis  and  exhibit  no  long-term 
stability.  This  class  of  compounds  (azides)  is  chemically  interesting  and  highly  energetic,  but  it 
exhibits  limited  structural  stability,  even  offering  explosive  character.  Further  research  in  this  area 
is  warranted,  however,  since  synthesis  techniques  for  aade-like  structures  are  often  applicable  to 
other  classes  of  compounds  that  exhibit  low  decomposition  barriers. 

4.3  Nitrated  Cyclic  and  Polycyclic  Compounds 

Nitramine  chemistry  is  presently  a  dynamic  area  for  synthesis  of  new  compounds.  Although 
compounds  such  as  RDX  and  HMX  have  been  known  for  decades,  there  was  little  new  activity 
until  the  1980s  when  interest  was  renewed  in  nitramines  based  on  the  development  of  caged 
polycyclic  structures.  These  new  compounds,  such  as  adamantane  and  wurtzitane,  are 
synthesized  by  forming  C-C  bridges  between  the  rings  of  two  or  more  alcane  structures.  The 
successful,  and  amazingly  simple,  synthesis  of  CL20  constituted  a  major  advance  in  synthesi^g 
nitramines  with  much  higher  energy  content.  Calculations  of  CL20  confirm  the  structure  of  this 
molecule  and  verify  the  high  positive  heat  of  formation.  Based  on  studies  of  CL20,  several  new 
polycyclic  nitramine  structures  were  analyzed.  The  first  compound:  3,5,8, 10, 1 1, 12-hexanitro- 
3,5,8,10,1 1,12-hexaazatetracyclo  □  (5.3.1.1.^’^0‘*’^)  dodecane,  is  formed  by  carbon-bridging  on 
two  RDX  molecules.  This  compound  exhibits  vibrational  stability  in  Cg  symmetry  but  has  a 
negative  oxygen  balance.  A  second  compound  of  this  type  can  be  formed  by  carbon-bridging  on 
two  HMX  molecules.  This  compound:  3,5,7,10,12,14,15,16-octanitro-3,5,7,10,12,14,15,16- 
octaazapentacyclo  □  (7.5.1.  l^’^O.'^’^^O.d,!!)  heaxadecane,  also  exhibits  vibrational  stability  in  Cg 
symmetry  but  still  has  an  oxygen  balance  of  -4  O  atoms.  Next,  nitrogen-bridging  on  two  RDX 
molecules  was  examined  and  a  stable  pentadecane  structure  was  found:  2,4,6,8,10,12,13,14,15- 
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nonanitra-2,4,6, 8,10,12, 13,14,15-nonaazatetracyclo  □  (5.5.1. pentadecane.  This  is  the 
first  polycyclic  nitramine  that  has  been  analyzed  that  exhibits  a  positive  oxygen  balance  (+3  0 
atoms).  Synthesis  activities  on  these  new  polycyclic  nitramines  appears  to  be  warranted. 

4.4  Dinitramide  Compounds 

As  is  now  well  known,  the  chemistry  of  dinitramide  compounds  was  developed 
independently  in  both  Russia  (Ref  27)  and  by  Bottaro  and  Schmitt*  at  SRI.  Several  new 
compounds,  mainly  alkali  and  ammonium  salts,  were  synthesized  that  exhibit  the  basic  N(N02)2 
structure.  UTRC  undertook  an  extensive  study  of  this  new  inorganic  anion  and  several  new 
dinitramide  compounds.  Calculations  confirmed  the  x-ray  structure  of  the  dinitramide  anion  and 
shed  light  on  decomposition  pathways.  Several  new  compounds  were  exaimned  for  structural 
stability,  including  N(N02)3  and  HADN.  Both  of  these  compounds  were  later  synthesized  by  SRI 
but  exhibit  limited  stability,  except  at  low  temperatures  [N(N02)3],  or  appear  to  hydrolyze  readily 
to  a  liquid  [HADN].  Extensive  propellant  formulation  activities  are  currently  in  progress  with  the 
parent  compound  of  this  group  [ADN],  since  it  is  a  proposed  chlorine-free  replacement  for  AP . 

This  dinitramide  group  can  also  be  realized  by  systematic  nitration  of  a  starting  compound 
such  as  hydrazine.  There  are  three  vibrationally  stable  compounds:  1,1-dinitrohydrazine,  1,2- 
dinitrohydrazine  and  tetranitrohydrazine.  The  first  and  last  compound  in  this  group  exhibit  a 
dinitramide  structure;  1,1-dinitrohydrazine  (aminodinitramide)  exhibits  the  -NH2  group, 
suggesting  that  the  activating  effect  of  this  primary  amine  group  could  lead  to  a  high  degree  of 
reactivity.  The  jewel  in  this  group  appears  to  be  tetranitrohydrazine,  which  exhibits  very  little 
central  N-N  bond  stretching  upon  complete  nitration  of  the  parent  hydrazine  molecule.  A 
performance  of  Igp  =  300  sec  is  predicted  using  this  oxidizer  \\dth  R45M  as  a  binder.  On  a  density 
impulse  basis,  nearly  15%  improvement  in  performance  can  be  realized.  This  compound  should 
be  a  leading  candidate  for  synthesis  activities  of  advanced  oxidizers. 

The  chemistry  of  the  nitroborane  [B(N02)2”]  structures  that  parallel  dinitramide  structures 
is  not  fully  understood.  Both  the  nitro  and  nitrite  forms  of  BH2NO2  are  conformationally  stable, 
but  there  appear  to  be  no  stable  boron  analogues  of  nitromethane  or  dinitramide.  Preliminary 
studies  indicate  that  B-bridged  structures  are  formed  but,  clearly,  more  studies  are  indicated  for 
this  chemistry  area. 

4.5  Nitrated  Ring  Structures 

Following  discussions  with  several  groups  working  on  the  development  of  explosives  with 
lessened  sensitivity  to  detonation,  UTRC  undertook  calculations  on  several  substituted  triazole 
and  tetrazole  compounds.  The  baseline  compound  that  was  examined  was  NTO  which  is  a 
prototype  as  an  insensitive  explosive.  Calculations  indicated  a  planar  structure  for  this 
compound  with  a  calculated  heat  of  formation  of -10  kcal/mol.  These  results  are  in  excellent 
agreement  with  the  known  experimental  characteristics  of  this  compound  and  suggest  that  this 
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level  of  theory  was  adequate  for  an  initial  evaluation  of  stability  and  energy  content.  The  several 
possible  decomposition  pathways  for  NTO  were  examined,  which  indicated  that  C-NO2  bond 
breaking  is  the  lowest  energy  decomposition  pathway.  Further  studies  of  the  decomposition 
pathways  for  NTO  are  currently  in  progress  by  K.  Lammertsma^  at  the  University  of  Alabama. 

Studies  were  then  focused  on  the  structure,  stability  and  energy  content  of  TNAZ.  There 
are  two  stable  conformations  for  this  compound  lying  only  0.3  kcal/mol  apart.  Further,  there  are 
very  low  barriers  for  nitro-group  rotation  in  this  system,  leading  to  a  complicated  spectra.  The 
calculated  heat  of  formation  for  TNAZ  is  +35.0  kcal/mol.,  a  value  much  higher  than  for  NTO. 
This  compound  should  be  further  studied,  including  an  accurate  analysis  of  its  thermochemistry. 

Beyond  TNAZ,  calculations  were  undertaken  on  several  compounds  that  were  suggested  or 
already  synthesized  by  SRI.  They  include  l-aniino-3-nitro- 1,2,4-triazole;  l-amino-3,5-dinitro- 
1, 2,4-triazole;  l-fluoro-3-nitro- 1,2,4-triazole;  and  1 -amino-3 -nitrotetrazole.  This  last  compound 
has  a  heat  of  formation  similar  to  CL20  [+96.5  kcal/mol]  and  exhibits  reasonable  stability.  The 
triazoles  that  were  examined  exhibit  +50  to  +70  kcal/mol  for  calculated  heats  of  formation. 

Finally,  l,3-dinitro-l,3-diazacyclobutane,  the  smallest  member  in  the  RDX,  HMX  nitramine 
series,  was  examined.  Studies  of  TNAZ  suggested  that  this  compound  should  exhibit  stability. 
This  compound  has  a  molar  heat  of  formation  similar  to  RDX  and  surprisingly  little  strain  energy 
in  the  ring.  HEDMX  is  suggested  as  an  identifier  for  this  new  cyclic  nitramine,  which  is  a  prime 
target  for  experimental  synthesis. 

4.6  Substituted  Prismanes 

The  last  series  of  compound,  [n]-prismanes,  was  not  completely  evaluated.  Experimental 
activities  at  Phillips  Laboratory  and  other  centers  indicate  that  this  class  of  compounds  may  be 
usefiil  in  improving  the  Igp  in  the  propellant  fuel/liquid  oxygen  (LP-l/LOX)  propulsion  systems. 
Cubane  addition  to  LP-1  has  already  been  demonstrated  to  yield  an  improvement  in  Ijp. 
Calculations  indicate  that  both  the  nitroprismanes  and  dinitroprismane  structures  are  stable. 
Nitration  lowers  the  heat  of  formation  by  ~  10%  but  these  compounds  exhibit  high  heats  of 
formation  (+100  to  +125  kcal/mol)  in  any  case,  owing  to  their  high  strain  energy.  Further 
synthesis  and  compound  characterization  activities  should  be  encouraged  for  this  class  of  strained 
hydrocarbons. 


^  Lanuneitsma,  K.:  private  communication.  University  of  Alabama,  1994. 
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INVESTIGATION  OF  TWO  POTENTIAL  VISIBLE  WAVELENGTH  CHEMICAL 
LASER  SCHEMES  BASED  ON  THERMAL  DISSOCIATION  OF  CHLORINE  AZIDE 

D.  J.  Benard,  B.  K.  Winker,  M.  A.  Chowdhury  and  T.  A.  Seder 
Rockwell  International  Science  Center 
Thousand  Oaks,  California  91360 

H.  H.  Michels 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 

ABSTRACT 

Chlorine  azide  and  fluorine  azide  were  synthesized  by  the  continuous  titration  of  dilute 
HN3  with  dilute  CIF  and  F2,  respectively.  The  reaction  products  were  analyzed  by  both  infrared 
and  ultraviolet  absorption  measurements  before  a  CO2  TEA  laser  was  used  to  pulse  heat  these 
energetic  molecules  in  the  presence  of  SFg  (sensitizer)  and  either  I2  or  CF3l/XeF2  gas  mixtures. 
The  resulting  thermal  dissociation  of  the  energetic  azides  produced  metastable  halogen  nitrenes 
which  reacted  by  energy  transfer  to  pump  either  the  NF(b)  or  IF(B)  states,  which  are  potentially 
capable  of  lasing  at  visible  wavelength.  Time  resolved  emission  spectroscopy  was  employed  to 
follow  these  emitters  as  well  as  the  metastable  intermediates  and  I*,  while  the  concentrations  of 
the  groimd  state  dihalides  were  tracked  by  laser  induced  fluorescence.  These  data  revealed  that 
the  global  rate  of  I*  quenching  was  negligible,  the  rate  of  NF(b)  quenching  by  NF(a,  X)  was  less 
than  6  x  10'^^  cm^/s,  the  yield  of  NCl(a)  per  CIN3  molecule  was  12%  and  the  radiative  rate  of 
the  NCl(a-X)  transition  was  measured  as  0.9/s  in  good  agreement  with  a  recent  theoretical 
estimate  by  Yarkony.  The  yield  of  NCl(a)  is  explained  by  ab  initio  calculations  of  the  CIN3 
potential  surfaces  which  demonstrate  that  a  singlet-triplet  crossing  occurs  inside  the  barrier  to 
dissociation.  Large  transient  concentrations  of  IF(X),  generated  (in  situ)  by  the  thermally 
initiated  reaction  of  XeF2  with  CF3I,  were  also  measured  absolutely  by  a  titration  technique  and 
(using  this  data)  the  effective  rate  of  pumping  of  the  interhalogen  to  its  B  state  (by  NF(a)/NCl(a) 
mixtures)  was  measured  as  6.5  x  10'^^  cm^/s.  This  low  rate  is  explained  by  NF(a)  pumping  out 
of  vibrationally  excited  levels  of  the  IF(X)  state  and  NCl(a)  pumping  into  vibrationally  excited 
levels  of  the  IF(B)  state  that  lie  above  the  predissociation  limit.  Thermal  dissociation  of  the 
stable  levels  of  the  IF(B)  state  is  also  likely  to  occur  at  rates  that  are  competitive  with  the 
spontaneous  radiation.  Additional  work,  focused  on  excitation  of  the  NS(B-X)  transition  by 
chemical  reactions  of  NF(a),  is  in  progress. 
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THEORETICAL  STUDIES  OF  SPIN-FORBIDDEN  RADIATIONLESS  DECAY  IN 
POLYATOMIC  SYSTEMS.  II.  RADIATIONLESS  DECAY  OF  a-N202. 
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Iowa  State  University 
Ames,  Iowa  50011 

J.  A.  Montgomery,  Jr.  and  H.  Harvey  Michels 
United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 

D.  R.  Yarkony 
The  Johns  Hopkins  University 
Baltimore,  Maryland  21218 

ABSTRACT 

The  stability  with  respect  to  spin-forbidden  radiationless  decay  of  the  previously  reported 
[J.  Chem.  Phys.  M,  7248  (1988)]  asymmetric  dimer  of  NO,  N-N-0-0  (a-N202)  is  considered. 
The  spin-allowed  decay  channel  a-N202  (^A')  —>  N20(X^S^  +  0(^D)  is  endoergic.  However, 
the  spin-forbidden  decay  channel  a-N202(^A')  N20(X^S'^  +  0(^P)  is  exoergic.  Large  scale 

multireference  configuration  interaction  wave  functions,  approximately  300,000  - 1,400,000 
configuration  state  functions,  based  on  double  zeta  polarization  and  triplet  zeta  polarization  bases 
are  used  to  study  this  process.  The  minimum  energy  crossing  of  the  ground  singlet  ^  A'  state  and 
the  lowest  excited  triplet  ^A"  state  was  determined  as  was  the  interstate  spin-orbit  coupling.  This 
electronic  structure  data  was  used  in  the  context  of  a  simple  one-  dimensional  model  to  show  that 
a-N202  is  rapidly  predissociated  to  N20(X^S''^  and  0(^P). 
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ON  THE  STRUCTURE  AND  THERMOCHEMISTRY  OF  HYDROGEN  DINITRAMIDE 

H.  H.  Michels  and  J.  A.  Montgomery,  Jr. 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 

ABSTRACT 

Five  stable  conformations  of  hydrogen  dinitramide  have  been  identified  through  ab  initio 
electronic  structure  calculations.  The  most  stable  has  a  secondary  amine  structure,  HN(N02)25 
and  a  calculated  heat  of  formation  (0  K)  of  28.4  kcal/mol.  Four  stable  inorganic  acid  structures 
are  found  to  lie  7.9  -  9.7  kcal/mol  higher  in  energy.  The  predicted  enthalpy  of  deprotonization  of 
N(N02)(NOOH)  is  AHacid(298  K)  =  302.6±2  kcal/mol,  making  this  one  of  the  strongest  known 
gas-phase  inorganic  acids.  We  estimate  the  electron  affinity  of  the  N(N02)2  radical  to  be 
5.0±0.5  eV  and  the  N-N  bond  dissociation  energy  of  HN(N02)2  to  be  48±8  kcal/mol.  It  is  shown 
that  unimolecular  decomposition  is  the  most  likely  low-energy  decomposition  mechanism  for 
hydrogen  dinitramide. 
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STRUCTURE  AND  STABILITY  OF  TRINITRAMIDE 

J.  A.  Montgomery,  Jr.  and  H.  H.  Michels 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 

ABSTRACT 

Ab  initio  calculations  of  the  structure  and  thermochemistry  of  trinitramide,  N(N02)3,  are 
reported.  A  vibrationally  stable  C3  structure  having  1 .545  A  N-N  bond  lengths  is  found  at  the 
MP2/6-3 IG*  level  of  theory.  The  heat  of  formation  of  trinitramide  is  estimated  to  be  59 
kcal/mol.  Thermal  decomposition  of  trinitramide  most  likely  occurs  via  N-N  bond  cleavage, 
which  is  estimated  to  require  26  kcal/mol. 
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ENERGY  TRANSFER  FROM  METASTABLE  NF  TO  BORON  HYDRIDE 

D.  J.  Benard  and  E.  Boehmer 
Rockwell  International  Science  Center 
Thousand  Oaks,  California  91358 

H.  H.  Michels  and  J.  A.  Montgomery,  Jr. 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 

ABSTRACT 

Gas  phase  reaction  of  B2H6  with  HN3  at  400°C  yields  tetraazidodiborane  (TADB),  which 
subsequently  dissociates  into  two  diazidoborane  molecules.  Exposure  of  TADB  to  vibrationally 
excited  SF6  generates  free  BH  radicals  that  may  be  pumped  by  energy  transfer  from  metastable 
NF(a^  A)  molecules,  obtained  by  simultaneous  dissociation  of  FN3.  Even  though  BH  radicals  are 
rapidly  consumed,  the  transient  yields  approach  2%  of  initial  donor  concentration.  The  rate 
constants  for  sequential  excitation  of  the  a^Il  and  A^H  states  of  BH  by  NF(a^A)  are  estimated  to 
be  ~8  X  10'^®  and  ~3  x  10'^  ^  cm^/s,  respectively,  and  transient  population  inversions  with  peak 
optical  gains  of  10%/cm  are  developed  on  the  BH(A-X)  transition  at  433  nm. 
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STRUCTURES,  BONDING,  AND  ENERGETICS  OF  N2O2  ISOMERS 

K.  A.  Nguyen  and  M.  S.  Gordon 
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J.  A.  Montgomery,  Jr.  and  H.  H.  Michels 
United  Technologies  Research  Center 
EastHartford,  Connecticut  06108 

ABSTRACT 


The  structures  and  energetics  of  the  N2O2  isomers  are  predicted  at  several  levels  of  theory. 
Both  single  reference  and  multireference  based  correlated  methods  were  used  to  determine  the 
structures  and  relative  energies.  Five  high-energy  minima  were  located  above  2NO  with 
QCISD(T)/6-3 1  l+G(2df)//MP2/6-3 1 1+G(d)  (PT2F/6-3 1  l+G(2df)//MCSCF/6-3  lG(d))  relative 
energies  of  ca.  38  (51),  46  (51),  61  (74),  69(74)  and  68  (80)  kcal/mol  for  l,2-diaza-3,4- 
dioxacyclobutene  (1),  bond  stretch  l,3-diaza-2,4-dioxa[1.1.0]bicyclobutane  (2),  a-N202  (3),  4, 
and  l,3-diaza-2,4-dioxa[1.1.0]bicyclobutane  (5),  respectively.  The  effect  of  basis  sets  on 
structures  is  small  within  a  given  level  of  theory.  The  MCSCF  structures  agree  reasonably  with 
those  of  MP2. 
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THEORETICAL  PREDICTION  OF  THE  STRUCTURES  AND  STABILITIES  OF 

AZIDAMINES 
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K.  O.  Christe 
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D.  A.  Dixon 
The  DuPont  Company 
Wilmington,  Delaware  19880 

ABSTRACT 

Ab  initio  molecular  orbital  and  density  functional  theory  calculations  show  that  azidamines, 
a  new  family  of  polynitrogen  compounds,  are  minima  on  their  respective  potential  energy 
surfaces.  The  geometries,  vibrational  frequencies,  and  heats  of  formation  are  predicted  for 
N(N3)3,  HN(N3)2,  the  N(N3)2'  anion,  and  the  N(N3)4‘^  cation.  All  of  these  compounds  are 
highly  energetic  materials  with  large  positive  heats  of  formation.  The  results  suggest  that  these 
polynitrogen  compounds  could  be  synthetically  accessible,  and  potential  methods  for  their 
syntheses  are  proposed. 
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fflGH  SPECIFIC  ENTHALPIES  FROM  THE  PHOTOCHEMICALLY  INDUCED 

ISOMERIZATION:  BOHoHBO 

J.  L.  Gole 
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H.  H.  Michels 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 

ABSTRACT 

The  photochemical  interconversion  of  BOH  to  the  considerably  more  stable  HBO  isomer  is 
evaluated.  The  BOH  and  HBO  isomers,  while  differing  considerably  in  energy,  possess  a 
significant  barrier  to  interconversion  on  the  ground  *  A'  potential  energy  surface.  However,  by 
accessing  the  low-lying  double  minimum  a  ^  A'  state,  we  define  an  interconversion  pathway  that 
facilitates  the  interconversion  process,  allowing  a  substantial  energy  release.  Here,  we  calculate 
the  potential  surfaces  for  the  X  *A'  and  a  ^A'  states  of  HBO  <-»  BOH  and  present  limited 
evaluations  for  the  b  ^A"  and  A  ^  A'  states.  We  outline  means  to  distinguish  the  BOH  and  HBO 
isomers,  following  the  progress  of  the  isomeri2ation  process  using  optical,  infrared,  and  mass 
spectrometric  techniques.  The  energetics  of  this  geometric  isomerization  process  and  the 
potential  utilization  of  BOH  as  a  storable  high  energy  material  are  discussed. 
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